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Abstract 

We summarize the standard factorization theorems for hard pro- 
cesses in QCD, and describe their proofs. 



1 Introduction 

In this chapter, we discuss the factorization theorems that enable one to apply 
perturbative calculations to many important processes involving hadrons. In 

*Originally published in A.H. Mueller, ed., "Perturbative QCD" , Adv. Ser. Direct. High 
Energy Phys. 5, 1 (1988). Authors' afSliations updated. 
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this introductory section we state briefly what the theorems are, and in 
Sects. EJtoIll we indicate how they are apphed in calculations. In subsequent 
sections, we present an outline of how the theorems are established, both in 
the simple but instructive case of scalar field theory and in the more complex 
and physically interesting case of quantum chromodynamics (QCD). 

The basic problem addressed by factorization theorems is how to calculate 
high energy cross sections. Order by order in a renormalizable perturbation 
series, any physical quantity is a function of three classes of variables with 
dimensions of mass. These are the kinematic energy scale(s) of the scatter- 
ing, Q, the masses, m, and a renormalization scale /i. We can make use of 
the asymptotic freedom of QCD by choosing the renormalization scale to be 
large, in which case the effective coupling constant g{fi) will be correspond- 
ingly small, g{fi) ~ 1 / ln{fi / Aqcb) ■ The renormalization scale, however, will 
appear in ratios Q/yU and yu/m, and at high energy at least one of these ratios 
is large. If we pick ~ Q, for instance, then at n loops the coupling will 
generally appear in the combination g'^'"'{Q) ln""((5/m), with a = 1 or 2. (See 
Sect. [7|) As a result, the perturbation series is no longer an expansion in a 
small parameter. The presence of logarithms involving the masses shows the 
importance of contributions from long distances, where the precise values of 
masses (including the vanishing gluon mass!) are relevant. For such contri- 
butions we do not expect asymptotic freedom to help, since it is a property 
of the coupling only at short distances. In summary, a general cross sec- 
tion is a combination of short- and long-distance behavior, and is hence not 
computable directly in perturbation theory for QCD. 

There are exceptions to this rule. For reasons which will become clear 
in Sect. [71 these are inclusive cross sections without hadrons in the initial 
state, such as the total cross section for e^e~ annihilation into hadrons, or 
into jets. 

This leaves over, however, the majority of experimentally studied lepton- 
hadron and hadron-hadron large momentum transfer cross sections, as well 
as inclusive cross sections in e^e~ annihilation with detected hadrons. Fac- 
torization theorems allow us to derive predictions for these cross sections, 
by separating (factorizing) long-distance from short- distance behavior in a 
systematic fashion. Thus almost all applications of perturbative QCD use 
factorization properties of some kind. 

In this chapter, we will explicitly treat factorization theorems for inclusive 
processes in which (1) all Lorentz invariants defining the process are large and 
comparable, except for particle masses, and (2) one counts all final states that 
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include the specified outgoing particles or jets. The second condition means 
that we consider such processes as hadron A + hadron B hadron C + X, 
where the X denotes "anything else" in addition to the specified hadron 
C. The first condition means that in this example the specified hadron C 
should have a transverse momentum comparable to the center-of-mass en- 
ergy. For such processes, the theorems show how to factorize long distance 
effects, which are not perturbatively calculable, into functions describing the 
distribution of partons in a hadron — or hadrons in a parton in the case 
of final-state hadrons. Not only can these functions be measured experi- 
mentally, but also the same parton distribution and decay functions will be 
observed in all such processes. The part of the cross section that remains 
after the parton distribution and decay functions have been factored out is 
the short distance cross section for the hard scattering of partons. This hard 
scattering cross section is perturbatively calculable, by a method which we 
describe below. 

Some examples of processes for which one expects a factorization theorem 
of this type to hold include (denoting hadrons hy A, B, C . . .) 

• Deeply inelastic scattering, lepton + A ^ lepton' -|- X; 

• e-^ + e- ^ A + X; 

• The Drell-Yan process, 

- A + B ^ + fj,- + X, 

- A + B^e+ + e-+X, 

- A + B + X, 

- A + B ^Z + X; 

• A + B ^jet + X; 

• A + B ^ heavy quark -|- X. 

In the last example, the heavy quark mass, which must be large compared 
to 1 GeV, plays the role of the large momentum transfer. In the Drell-Yan 
case, the kinematic invariants are the particle masses, the square, s, of the 
center-of-mass energy, and the invariant mass Q and transverse momentum 
gj_ of the lepton pair. The requirement, for the theorems that we discuss, 
that the invariants all be large and comparable means that not only should 
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be of order s, but also that either we integrate over all or qx is of order 

Q. 

There are applications of QCD to processes in which there is a large 
momentum scale involved but for which the most straightforward sort of fac- 
torization theorem, as discussed in this chapter, must be modified. However, 
the same style of analysis as we will describe applies to these more general sit- 
uations. (The Drell-Yan process when qx is much less than Q is an example.) 
We will summarize these in Sect. ITUl 

Some of the factorization properties, such as those we describe in this 
chapter, have been proved at a reasonable level of rigor within the context of 
perturbation theory. But many of the other results have, so far, been proved 
less completely. 

The following three subsections give explicit factorization theorems for 
three basic cross sections from the list above, deeply inelastic scattering, 
single-particle inclusive annihilation and the Drell-Yan process. These three 
examples illustrate most of the issues involved in the application and proof 
of factorization. We close the section by relating factorization to the parton 
model. 



1.1 Deeply Inelastic Scattering 

Deeply inelastic lepton scattering plays a central role in any discussion of fac- 
torization, both because this was the first process in which pointlike partons 
were "seen" inside the hadron, and because much of the data that determines 
the parton distribution functions comes from measurements in this process. 
In particular, let us consider the process e + A ^ e + X , which proceeds via 
the exchange of a virtual photon with momentum q^. From the measured 
cross section, one can extract the standard hadronic tensor W^^{q^,p^), 

W^'' = ^ I dSe'''-yY.{A\f{y)\X){X\3''mA) 

+P2\X,(^ ) , [i) 

p-q 

where = —q^q^, x = Q'^/2q-p, p^ is the momentum of the incoming hadron 
A, and j'^(x) is the electromagnetic current. (More generally, j'^iy) can be 
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any electroweak current, and there will be more than two scalar structure 
functions Fj.) 

We consider the process in the Bjorken limit, i.e., large Q at fixed x. The 
factorization theorem is contained in the following expression for W^^ , 

W^''{q^,pn = Y. /'^/aM(e,/^)^^a^V,eP^/^,«s(/x))+ remainder. (2) 

Here fa/A{,i^lj) is a parton distribution function, whose precise definition is 
given in Sect.lH There, /a/A(^,Ai)d^ is interpreted as the probability to find 
a parton of type a (= gluon, u, u, d, d, . . . ) in a hadron of type A carrying a 
fraction ^ to ^ + d^ of the hadron's momentum. In the formula, one sums 
over all the possible types of parton, a. We can prove Eq. (j21) in perturbation 
theory, with a remainder down by a power of Q (in this case, the power is 
Q~'^ modulo logarithms, but the precise value depends on the cross section 
at hand, and has not always been determined). 

We can project Eq. Q onto individual structure functions: 

-J- fa/A{^,l^) Hia 1^-, — + remainder, 
-F2{x,Q^) = ^ / ^ /a/A(^,/i) -H2a ( ) + remainder, (3) 

The extra factors of 1/x and ^/x in the equation for F2 are needed because 
of the dependence on target momentum of the tensor multiplying F2. 

Inspired by the terminology of the operator product expansion for the 
moments of the structure functions, it is conventional to call the first term 
on the right of either of Eqs. or the leading twist contribution, and 
to call the remainder the higher twist contribution. The same terminology 
of leading and higher twist is used for the factorization theorems for other 
processes. 

It is not so obvious why proving Eq. (j2)) in perturbation theory is useful, 
given that hadrons are not perturbative objects. But suppose we do decide 
on a way of computing the matrix elements in Eq. perturbatively. For any 
such formulation for hadron A, both W'^'^ and fa/A will depend on phenomena 
at the scale of hadronic masses (or some other infrared cutoff), and the exact 
nature of these phenomena will depend on our particular choice of A, as well 
as on the precise values we pick for both hadronic and partonic masses. The 
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content of the factorization theorem is that this dependence of W'^'^ on low 
mass phenomena is entirely contained in the factor of fa/A- 

The remaining function, the hard scattering coefficient Hj^'^, has two im- 
portant properties. First, it depends only on the parton type a, and not 
directly on our choice of hadron A. Secondly, it is ultraviolet dominated, 
that is, it receives important contributions only from momenta of order Q. 
The first property allows us to calculate Hj^'^ from Eq. Q with the simplest 
choice of external hadron, A = b, b being a parton. (We will see an example 
of this in our calculations for the Drell-Yan process in Sect. El) The second 
property ensures that when we do this calculation, Hl^^ will be a power se- 
ries in as{Q), with finite coefficients. We now assume that nonperturbative 
long-distance effects in the complete theory factorize in the same way as do 
perturbative long-distance effects. Once this assumption is made, we can 
interpret our perturbative calculation of Hj^'^ as a prediction of the theory. 
Parton model ideas, summarized in Sect 11.41 give motivation that the as- 
sumption is valid. Note that our definition of the parton distributions, which 
we will give in Sect. |31 is an operator definition, which can be applied beyond 
perturbation theory. 

This ability to calculate the Hj^'^ results in great predictive power for fac- 
torization theorems. For instance, if we measure F2{x,Q^) for a particular 
hadron A, Eq. Q will enable us to determine fa/ a- We then derive a predic- 
tion Q^) for the same hadron A, in terms of the observed F2 and the 
calculated functions Hia- This is the simplest example of the universality of 
parton distributions. 

The functions Hia may be thought of as hard-scattering structure func- 
tions for parton targets, but this interpretation should not be taken too 
literally. In any case, methods for putting this procedure into practice, in- 
cluding definitions for the parton distributions are the subjects of Sects. |21to 



Originally, Eq. Q was primarily discussed in terms of the moments of 
the structure functions, such as 
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With this notation, Eq. becomes 

Fi{n, Q^) = ^ fa/Ain, fi) hJti, — , ^^(/i) J . (5) 

In this form of the factorization theorem, when n is an integer, the fj/A{n, /i) 
are hadron matrix elements of certain local operators, evaluated at a renor- 
malization scale fi. On the other hand, the structure function moment 
F2{n, Q"^) can be expressed in terms of the hadron matrix element of a prod- 
uct of two electromagnetic current operators evaluated at two nearby space- 
time points. Equation (0) thus appears as an application of the operator 
product expansion [11121 El- The product of the two operators is expressed 
in terms of local operators and some perturbatively calculable coefficients 
Hia{n,Q/ fi,as{^)), called Wilson coefficients. It was using this scheme that 
the Hia{n,Q/ fi,as{fi)) were first calculated 



1.2 Single Particle Inclusive Annihilation 

In this subsection, we consider the process 7* —>■ A + X , where 7* is an 
off-shell photon. The relevant tensor for the process, for which structure 
functions analogous to those in Eq. (P) may be derived, is 

D^''{x,Q) = ^ld'y e^^-^5^(0|j^(i/)|AX)(AX|j-(0)|0), (6) 

where is now a time-like momentum and = q^. The sum is over all 
final-states that contain a particle A of defined momentum and type. We 
define a scaling variable hy z = 2p ■ qjQ^ ^ where j)^ is the momentum of A, 
and we will consider the appropriate generalization of the Bjorken limit, that 
is, Q large with z fixed. 

The factorization theorem here is quite analogous to Eq. (|21), but incor- 
porates the slightly different kinematics, 

D^\z, Q) = J2 ['t ^-^'(^/^' ^/^^ "^(/^)) ^^/-^(O, (7) 

with corrections down by a power of Q, as usual. We have used the same 
notation for the hard functions as in deeply inelastic scattering, and as in that 
case they are perturbatively calculable functions. Here it is the fragmentation 
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functions dA/a{0 which are observed from experiment, and which occur in 
any similar inclusive cross section with a particular observed hadron in the 
final state. For example, single-particle inclusive cross sections in deeply 
inelastic scattering cross sections require the factorization both of parton 
distributions fa/A^ with A the initial hadron, and of distributions ds/a^ with 
B the observed hadron in the final state. We shall not go into the details of 
such cross sections here jS]. 



1.3 Drell-Yan 

Our final example to illustrate the important issues of factorization is the 
Drell-Yan process: 

A + 5 ^ /i+ + /i- + X (8) 

at lowest order in quantum electrodynamics but, in principle, at any order in 
quantum chromodynamics. is now the momentum of the muon pair. We 
shall be concerned with the cross section dcr/dQ^d?/, where is the square 
of the muon pair mass, 

= q'q,, (9) 
and y is the rapidity of the muon pair, 

y = h-('-^)- (10) 



2 \1-Pb 

We imagine letting and the center of mass energy ^/s become very large, 
while Q'^/s remains fixed. 

The relevant factorization theorem, accurate up to corrections suppressed 
by a power of Q^, is 



da 



_ I, J XA J XB 

X fa/A{U,t^) Habl^,^,Q;y:,as{^J')) fb/B{^B,^J')- (H) 

V4a 4b W / 



Here a and b label parton types and we denote 



'02 /r)2 
a;^ = e^/^, XB = e-yJ^. (12) 
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The function Hab is the ultraviolet-dominated hard scattering cross sec- 
tion, computable in perturbation theory. It plays the role of a parton level 
cross section and is often written as 

da 

when it is not necessary to display the functional dependence of Hab on the 
kinematical variables. The parton distribution functions, /, are the same as 
in deeply inelastic scattering. Thus, for instance, one can measure the parton 
distribution functions in deeply inelastic scattering experiments and apply 
them to predict the Drell-Yan cross section. As before, the parameter is a 
renormalization scale used in the calculation of Hab- 



1.4 Factorization in the Parton Model 

Having introduced the basic factorization theorems, we will now try to give 
them an intuitive basis. Here we shall appeal to Feynman's parton modellHj. 
In fact, we shall see that factorization theorems may be thought of as field 
theoretic realizations of the parton model. 

In the parton model, we imagine hadrons as extended objects, made up of 
constituents (partons) held together by their mutual interactions. Of course, 
these partons will be quarks and gluons in the real world, as described by 
QCD, but we do not use this fact yet. At the level of the parton model, 
we assume that the hadrons can be described in terms of virtual partonic 
states, but that we are not in a position to calculate the structure of these 
states. On the other hand, we assume that we do know how to compute the 
scattering of a free parton by, say, an electron. By "free", we simply mean 
that we neglect parton-parton interactions. This dichotomy of ignorance 
and knowledge corresponds to our inability to compute perturbatively at 
long distances in QCD, while having asymptotic freedom at short distances. 

To be specific, consider inclusive electron-hadron scattering by virtual 
photon exchange at high energy and momentum transfer. Consider how this 
scattering looks in the center-of-mass frame, where two important things 
happen to the hadron. It is Lorentz contracted in the direction of the colli- 
sion, and its internal interactions are time dilated. So, as the center-of-mass 
energy increases the lifetime of any virtual partonic state is lengthened, while 
the time it takes the electron to traverse the hadron is shortened. When the 
latter is much shorter than the former the hadron will be in a single virtual 
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state characterized by a definite number of partons during tlie entire time tlie 
electron takes to cross it. Since tlie partons do not interact during this time, 
each one may be thought of as carrying a definite fraction x of the hadron's 
momentum in the center of mass frame. We expect x to satisfy < x < 1, 
since otherwise one or more partons would have to move in the opposite di- 
rection to the hadron, an unlikely configuration. It now makes sense to talk 
about the electron interacting with partons of definite momentum, rather 
than with the hadron as a whole. In addition, when the momentum trans- 
fer is very high, the virtual photon which mediates electron-parton scattering 
cannot travel far. Then, if the density of partons is not too high, the electron 
will be able to interact with only a single parton. Also, interactions which 
occur in the final state, after the hard scattering, are assumed to occur on 
time scales too long to interfere with it. 

With these assumptions, the high energy scattering process becomes es- 
sentially classical and incoherent. That is, the interactions of the partons 
among themselves, which occur at time-dilated time scales before or after the 
hard scattering, cannot interfere with the interaction of a parton with the 
electron. The cross section for hadron scattering may thus be computed by 
combining probabilities, rather than amplitudes. We define a parton distri- 
bution fa/niO the probability that the electron will encounter a "frozen", 
noninteracting parton of species a with fraction ^ of the hadron's momentum. 
We take the cross section for the electron to scatter from such a parton with 
momentum transfer as the Born cross section (Tb{Q'^,0- Straightforward 
kinematics shows that for free partons ^ > x = 2p- q/Q'^, and the total cross 
section for deeply inelastic scattering of a hadron by an electron is 

aeH{x,Q') = fdUa/H{OM^/^,Q')- (14) 

This is the parton model cross section for deeply inelastic scattering. It is 
precisely of the form of Eq. (j2I), and is the model for all the factorization 
theorems which we discuss in this chapter. 

Essentially the same reasoning may be applied to single-particle-inclusive 
cross sections and to the Drell-Yan cross section. For example, in the par- 
ton model the latter process is given by the direct annihilation of a parton 
and anti-parton pair, one from each hadron, in the Born approximation, 
ct'qIQ'^ ,y). The interactions which produce the distributions of each such 
parton occur on a scale which is again much longer than the time scale of the 
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annihilation and, in addition, final-state interactions between the remaining 
partons take place too late to affect the annihilation. We thus generalize (|14jl 
to the parton model Drell-Yan cross section 



da 



f f d^B fa/A{U) f-a/B{iB) a'j,{Q\y), (15) 

„ JxA JxR 



where xa,b are defined in^J Equation (fT3j) is of the same form as the full 
factorization formula, (|TT|l . except that there is only a single sum over par- 
ton species, since the hard process here consists of a simple quark-antiquark 
annihilation. In the parton model, the functions fa/A{^A) in Drell Yan must 
be the same as in deeply inelastic scattering, Eq. (HH), since they describe 
the internal structure of the hadron, which has been decoupled kinemati- 
cally from the annihilation and from the other hadron. It is important to 
notice that the Lorentz contraction of the hadrons in the center of mass sys- 
tem is indispensable for this universality of parton distributions. Without 
it the partons from different hadrons would overlap a finite time before the 
scattering, and initial-state interactions would then modify the distributions. 

We now turn to the technical discussion of factorization theorems in QCD, 
but it is important not to loose sight of their intuitive basis in the kinematics 
of high energy scattering. In fact, when we return to proofs of factorization 
theorems in gauge theories (Sects. |H1 and ^ these considerations will play a 
central role. 



2 Calculation of the Hard Scattering Cross 
Section 

In this and the following two sections, we discuss the explicit calculation of 
the hard scattering functions for the Drell Yan cross section. In doing so, 
we will cover most of the technical points which are encountered in applying 
factorization in other realistic cases as well. 

At order zero in as for the Drell-Yan cross section, the hard process 
described by Hab is quark-antiquark annihilation, as illustrated in^J One can 
simply compute this parton level cross section from the Feynman diagram 
and insert it into the factorization formula (jllj) . The resulting cross section 
is not itself a prediction of QCD, although it is a prediction of the parton 
model. The factorization theorem will make the connection between the two. 
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Figure 1: Bom amplitude for the Drell-Yan process. 
At the Born level, it is natural to define /a/a(0 = ^{^ ~ 0- We then find 

where the factor 5^1 indicates that parton a must be the antiparticle to 
parton b. Here C(e) is 1 if we work in 4 space-time dimensions. However, 
when one wants to calculate higher order contributions, it will turn out to 
be useful to perform the entire calculation in 4 — 2e dimensions. Then 

W J W J (l-2e/3)(l-2e)r(l-2e)- 

The e dependence here arises from three sources. First, the Dirac trace 
algebra gives an angular dependence 1 + cos^ 9 — 2e. Secondly, one introduces 
a factor {fj,'^ / {Att) e'^Y so as to keep the cross section at a constant overall 
dimensionality^ of M~'^. Finally, the integration over the lepton angles in 
4 — 2e dimensions gives the remaining e dependence. Actually, it is quite 
permissible to perform the lepton trace calculation and the integration over 
lepton angles in 4 dimensions instead of 4 — 2e dimensions. This procedure 
results in multiplying the Born cross section and the higher order cross section 
by a common, e-dependent factor. As we will see below, such a factor will 
drop out in the physical cross section. 

Now let us calculate H at one loop. At first order in a^, the cross 
section gets contributions from the graphs shown in |21 along with their 

^We use (/i^/(47r) e^)^ rather than (/i^)^ in anticipation of our use of MS renormaliza- 
tion. 
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Figure 2: Order contributions to the Drell-Yan cross section. 

mirror diagrams. In this figure, we show contributions to both the am- 
phtude and its complex conjugate, separated by a vertical line which rep- 
resents the final state. We will use this notation frequently below, and 
refer to diagrams of this sort as "cut diagrams". The situation now is 
not so simple, because a straightforward calculation of the cross section for 
quark + antiquark ^ fi'^ -\- fi" + X according to the diagrams shown above 
yields an infinite result when we use massless, on-shell quarks as the incoming 
particles. 

Following Sect. 11.11 we use the factorization formula (lllj) applied to in- 
coming partons instead of incoming hadrons. Since the details associated 
with parton masses are going to factorize, we can choose to calculate the 
cross section for parton a + parton b ^ fi~^ + fi~ +X with the partons having 
zero mass and transverse momentum. Let us call this cross section Gab'- 



In this calculation there are both ultraviolet and infrared divergences. Di- 
mensional regularization is used to regulate them both. The factorization 
formula is then 




(18) 




(19) 
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Both factors in the formula depend on /i, which is the scale factor introduced 
in the dimensional regularization and subsequent MS renormalization |lj of 
Green functions of ultraviolet divergent operators. One introduces a factor 

for each integration J d'^~'^'^k in order to keep the dimensionality of the result 
independent of e. Ultraviolet divergences then appear as poles in the variable 
e, which are subtracted away, as explained in [7j. The factor e'^/^in) that 
comes along with the fi is the difference between MS renormalization and 
minimal subtraction (MS) renormalization. Here 7 = 0.577... is Euler's 
constant. 

Let us suppose that we have calculated Gab to two orders in perturbation 
theory. We denote the perturbative coefficients by 

Gab = G^; + ^G^^^ + Oial). (21) 

TT 

Thus G^^fj is the Born cross section in Eq. p6|). G^^ij is the first correction. 

The first correction G^^ij will generally have ultraviolet divergences at 
e = 0, coming from virtual graphs, and these divergences will appear as 
1/e poles. Following the minimal subtraction prescription, we remove these 
ultraviolet poles as necessary.^ In general 1/e poles of infrared origin will 
remain in and we shall discuss these infrared poles presently. 

Let us similarly denote the perturbative coefficients of the hard scattering 
cross section Hat by 

^.. = ^i? + -^i;^ + (^(«D- (22) 

71 

It is these coefficients that we would like to calculate. 

All we need to know to calculate H from G is the perturbative expansion 
of the functions fa/b{x,e), which, according to the factorization theorem, 

^In the particular case of the Drell-Yan cross section (or, more generally, a cross section 
for which the Born graph represents an electroweak interaction), the first QCD correction 

^ab '^''^ ultraviolet divergent, provided that we include the propagator correc- 

tions for the incoming quark lines This follows from (1) the Ward identity expressing the 
conservation of the electromagnetic current and (2) the fact that the photon propagator 
does not get strong interaction corrections, at lowest order in QED. It can also be verified 
easily by explicit computation. 
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contain all of the sensitivity to small momenta, and are interpreted as the 
distribution of parton a in parton h. These functions can be calculated in a 
simple fashion using their definitions (Sect. |3]) as matrix elements (here in 
parton states) of certain operators. When the ultraviolet divergences of the 
operators are also renormalized using minimal subtraction, one finds simply 

/a/.(x; e) = 5,, ^(1 - ^) - ^ f + 0{al), (23) 

where P|^^(a;) is the lowest order Altarelli-Parisi [HJ kernel that gives the 
evolution with /i of the parton distribution functions. We will discuss the 
computations that lead to Eq. in Sect. IH For now, let us assume the 
result. 

When we insert these perturbative expansions ()23|) into the factorization 
formula, we obtain 

G^ab (^A, Xb, Q; ^' ^) + ~ ^ab (^^A, Xb, Q] 

= ^1? (^^A, xb, Q; ^; + ^ H^l^ (^xa, xb, Q; ^; 

(24) 

We can now solve for Hat- At the Born level, we find 

H^^ (^XA, XB, Q; ^; = (^xa, xb, Q; |; ■ (25) 
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Then at the one loop level we obtain 

Gab (^xa,xb,Q; 

(26) 

Thus the prescription is quite simple. One should calculate the cross sec- 
tion at the parton level, G^^/j , and subtract from it certain terms consisting of 
a divergent factor 1/e, the Altarelli-Parisi kernel, and the Born cross section 
(with e 7^ 0). The result is guaranteed to be finite as e — 0. 

Recall that the Born cross section G^^^ consists of an e dependent factor 
C(e) times the Born cross section in 4 dimensions, where C(e) arises from such 
sources as the integration over the lepton angles in the Drell-Yan process. A 
convenient way to manage the calculation is to factor C(e) out of the first 
order cross section G^^^ also. Then the prescription is to remove the 1/e pole 
in G^^\e)/C{e), set e = 0, and multiply by C(0) = 1. Thus we see that a 
function of e that is a common factor to and G^^^ cancels in the physical 
hard scattering cross section, as was claimed after Eq. (fTTj) . 

When calculating G^^\ it should be noted that there are contributions 
involving self energy graphs on the external lines, as in |21 The total of all 
external line corrections gives a factor of ^/z2 for each external quark (or 
antiquark) line and y/z^ for each external gluon line. Here Z2 and Z3 are the 
residues of the poles in the renormalized quark and gluon propagators. In 
the massless theory these have infrared divergences. For example the value 
of Z2 in massless QCD in Feynman gauge is 

^2 = l + ^ + 0(«D- (27) 

Then the contribution of the self energy graphs to G^^^ is a factor 2as/37re 
times the Born cross section. 



H 



(1) 

ab 



Xa 
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3 Relation to the renormalization group 



The prescription ()26|) for removing infrared poles is intimately related to the 
/X dependence of H^^^ — that is, to the behavior of if^^'* under the renor- 
malization group. In this section, we display this connection and show how 
it leads to the approximate invariance of the computed cross section under 
changes of (Of course, the complete cross section, to all orders of per- 
turbation theory, is exactly invariant under changes of yU. What we are now 
concerned with is the behavior of a finite-order approximation.) 

We recall that the Born cross section = H^^j contains some ^ depen- 
dence from the factor C{fi/Q,e), as specified in Eq. flT^ . The one loop cross 
section G^^^ contains this same factor, and we can simply factor it out of Eq. 
(j2ni) and set it to 1 when we set e = at the end. In addition, G^^jj contains 
a factor /i^^ from the loop integration. 

The {e^ jjL^ /Any factor multiplies the 1/e poles in G^l''- Writing 

j/i'^ = ^ + 2A ln(/i) + 0(e), (29) 

and reading off the value of A from Eq. (j26|l . we find the /i dependence of 
Gil^-and thus oiH^^: 



Q 



id i ^A, , Qj ■ 

(30) 



Here we have set e = and have suppressed the notation indicating e depen- 
dence; we have also noted that H^^^ does not depend on n when e = 0, so we 
have suppressed the notation indicating fi dependence in 
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We see that H^^^ contains logarithms oi fi/Q. If /i is fixed while Q becomes 
very large, then these logarithms spoil the usefulness of perturbation theory, 
since the large logarithms can cancel the small coupling that multiplies 
H^^\ For this reason, one chooses /i such that ln(yu/Q) is not large. For 
example, one chooses = Q or perhaps fi = 2Q or = Q/2. 

The freedom to choose /i results from the renormalization group equations 
obeyed by H and fa/A^O- The renormalization group equation for Hab is 



1^-^ Hab (^xa, Xb, Q; ^, 



~^J^ (^CAPc/a{CA,asifx)) Hcb(^^,XB,Q;-^,as{fx)^ 
- / dCs Pd/biCs, asifJ')) Had { XA, 7^, Q 

^ JxB V 



(31) 



1 ^^"^(/i) ) • 



Here Pc/a(C, «s(/^)) is the all orders Altarelli-Parisi kernel. It has a pertur- 
bative expansion 



Pc/a(e,«.(/i)) = ^pi;i(0 + --- (32) 



vr 

5(1) 



where -P!'/a(0 is the function that appears in Eq. (j2SI)- Thus at lowest order 
the renormalization group equation (jHTj) is a simple consequence of differen- 
tiating Eq. 

Parton distribution functions also have a dependence, which arises from 
the renormalization of the ultraviolet divergences in the products of quark 
and gluon operators in the definitions of these functions, given in Eqs. (I43|) 
and (j44|) below. The renormalization group equation for the distribution 
functions is 

^Jil ^"'^^^'^^ ^ Efc/g^ f ^a/fe(C,«s(/i)) /6M(|,/i) • (33) 

The physical cross section does not, of course, depend on fi, since fi is 
not one of the parameters of the Lagrangian, but is rather an artifact of 
the calculation. Nevertheless, the cross section calculated at a finite order 
of perturbation theory will acquire some /i dependence arising from the ap- 
proximation of throwing away higher order contributions. To see how this 
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comes about, we differentiate Eq. (fTT|) with respect to /i and use Eqs. PT|1 
and (jnS))- This gives 



d da 



dfi dQ'^dy 



a,b,c 



U (a Jx 



'XA -^KA 



1 rl rl 

J 1* 



fa/A {U, /^) ^c/a(CA, ".(/i)) H^b i — 7^, Q; 7:, I fb/si^B, 

+ i? terms. 

(34) 

Here the two terms shown relate to the evolution of the partons in hadron 
A. As indicated, two similar terms relate to the evolution of the partons in 
hadron B. We now change the order of integration in the second term to 
put the integration inside the Ca integration, then change the integration 
variable from to = ^a^a-, and finally reverse the order of integrations 
again. This gives 

d da 



d/x dQ^dy 

abc ''^^ -^^A ^A J^B 



-'Ca Jx 



B 



a,b,c 



X •^'^M^^' ^c/a(CA, as(Ai)) If' ^' ^' fb/B{^B, /i) 

+ 5 terms. 

(35) 
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We see that the two terms cancel exactly as long as Pa/b and Hab obey the 
renormalization group equations exactly. Now, when Hab is calculated only 
to order , it only obeys the renormalization group equation (jHT|) to the 
same order. In this case, we will have 

when the parton distribution functions obey the renormalization group equa- 
tion with the Altarelli-Parisi kernel calculated to order or better. One 
thus finds that the result of a Born level calculation can be strongly /i de- 
pendent, but by including the next order the /x dependence is reduced. 

We have argued that one should choose /i to be on the order of the large 
momentum scale in the problem, which is Q in the case of the Drell-Yan cross 
section. We have the right to choose n as we wish because the result would 
be independent of /i if the calculation were done exactly. The choice /i ~ Q 
eliminates the potentially large logarithms in Eq. ()30p . Another choice is 
often used. One substitutes for in Eq. (fTTj) the value = V^aCbS- 
We now have a value of yU that depends on the integration variables in the 
factorization. 

Let us examine whether this is valid, assuming that P^/b and Hab are 
calculated exactly. We replace /z by 

^l{\U^B)=^il~^ (^/UiBs)\ o<A<i. (37) 

At A = we have a valid starting point. When we get to A = 1 we have 
the desired ending point. The question is whether the derivative of the cross 
section with respect to A is zero. Applying the same calculation as before, 
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we obtain instead of Eq. the result 

ax doM? = £ L cT 1 2 j 

X ^a/c(CA, a.(/i(A, ^A, -^b))) /c/A /i(A, ^A, '^b)^ 

X /a/A (^A./ilAj^A.-^B)) Pc/a (Ca, as(/i(A, ^A, '^s))) 

X Ed, ( rr-, 7^, Q; I fb/si^B, /i(A, ^A, ^s)) 

V4aCa 4b V / 

+ i? terms. 

(38) 

Now making the same change of variables as before, we obtain 
d da 



dX dQ^dy 

j2 rdar^rd«4>"^'^'" 



X Pa/c{.CA,(^s{.IJ'{.\U,iB))) fc/A{^,IJ'{.\U,iB)^ 

^ rr (XA Xb ^ /i(A,^A,^B) ( ^ f (c I \ C C \\ 

\4a 4b V / 

X fa/A (^,K\U/Ca,^b)^ Pc/aiCA,as{fx{\,U/CA,^B))) 

X i^cfe 7^ ,a.(/i) fb/B{^B,l^{KU/CA,^B)) 

\4a 4b W / 

+ B terms. 

(39) 
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We see that the cancellation between the two terms has been spoiled, first by 
the differences in the values of /i(A, . . . ) in the two terms, but more impor- 
tantly by the differences in the arguments of the logarithm in the two terms. 
We conclude that the substitution of s for /x^ results in an error of order 
as no matter how accurately the hard scattering cross section is calculated. 
This is not a problem if the hard scattering cross section is calculated only 
at the Born level, which is, in fact, commonly the case. However, it is wrong 
to substitute s for fi"^ when a calculation beyond the Born level is used. 

4 The parton distribution functions 

The parton distribution functions are indispensable ingredients in the fac- 
torization formula We need to know the distribution of partons in a 
hadron, based on experimental data, in order to obtain predictions from the 
formula. In addition, we need to know the distribution of partons in a parton 
in order to calculate the hard scattering cross section Hab- The hard scatter- 
ing cross section is obtained by factoring the parton distribution functions 
out of the physical cross section. Evidently, the result depends on exactly 
what it is that one factors out. 

4.1 Operator Definitions 

In this section, we describe the definition for the parton distribution functions 
that we use elsewhere in this chapter. A more complete discussion can be 
found in Ref. In this definition, the distribution functions are matrix 
elements in a hadron state of certain operators that act to count the number 
of quarks or gluons carrying a fraction ^ of the hadron's momentum. We state 
the definition in a reference frame in which the hadron carries momentum 

with a plus component P~^, a minus component = m^/2P+, and 
transverse components equal to zero. (We use P± = (P° ± P^)/V2). 

The definition may be motivated by looking at the theory quantized on 
the plane x"*" = in the light-cone gauge A'^ = 0, since it is in this picture 
that field theory has its closest connection with the parton model 1Q\ . In this 
gauge, Q = 1, where ^ is a path-ordered exponential of the gluon field that 
appears in the definition of the parton distributions. The light-cone gauge 
tends to be rather pathological if one goes beyond low order perturbation 
theory, and covariant gauges are preferred for a complete treatment. How- 
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ever quantization on a null plane in the light-cone gauge provides a useful 
motivation for the complete treatment. 

In this approach the quark field has two components that represent the 
independent degrees of freedom; 'y~^ip{x) contains these components and not 
the other two. One can expand the two independent components in terms 
of quark destruction operators b{k~^, s) and antiquark creation operators 
d{k~^, k_i, sY as follows: 



+ , / X 1 ^ r dk+ r 



dki 



X [^+U{k, s)e-'''-^b{k+ , k^, s) + -f+V{k, s)e+'''-''d{k+, k±, s)^] . (40) 

The quark distribution function is just the hadron matrix element of the 
operator that counts the number of quarks. 

/.m(0 = ^ E / dk^ {P\ b{CP\ k^, s)H^P\ k^, s) IP). 

(41) 

In terms of ip{x), this is 

fMO = ^ j dx-e-^«^^^'"(P|V^(0,a;-, 0^)7+^(0, 0,Ox)|P). 

We can keep this same definition, while allowing the possibility of computing 
in another gauge, by inserting the operator 

g = Ve^V d2/-A+(0,i/-,0x)tc| , (42) 

where V denotes an instruction to order the gluon field operators A^{^, y^, 0_l) 
along the path. The operator Q is evidently 1 in the = gauge. With 
this operator, the definition is gauge invariant. 
We thus arrive at the definition [HI E] 

MO = ^J dx-e-^«^^^'(P|^(0,x-,0^)7+^V^(0,0,0x)|P). (43) 
For gluons, the definition based on the same physical motivation is 

/,m(0 = ^^ / dx-e-*«^^^-(P|P,(0,x-,0x)+''^a.F6(0,0,0^)/|P), 

(44) 
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where is the gluon field strength operator and where in Q we now use 
the octet representation of the SU(3) generating matrices tc- 



4.2 Feynman rules and eikonal lines 

The Feynman rules for parton distributions are derived in a straightforward 
manner from the standard Feynman rules. Consider, for instance the distri- 
bution fg/g of a quark in a quark. To compute this quantity in perturbation 
theory, we use the following identity satisfied by any ordered exponential, 

V expi^ig dXn ■ A{{X + r])n'')^ ^ P expi^ig dAn ■ v4(An'')| . 



(45) 



Using in Eq. 
states and write 



for instance, enables us to insert a complete set of 

M0 = ^ I dx-e-«^^^" J](P|f^(0,a:-,0x)|n)7+(n|vi/(0,0,0x)|P) , 

(46) 

where we define \E' as the quark field times an associated ordered exponential. 



\E'(x) = ipi^x) V exp 



dXv- A{x + Xv^) 



(47) 



where = gt-, and A^{x) = A'^{x)tc- To express the matrix elements in 
Eq. pUj) in terms of diagrams, we note that by (P7j) the gluon fields in the 
expansion of \l/ are time ordered by construction. Expanding the ordered 
exponentials, and expressing them in momentum space we find 



V exp 



dXn-A{Xn'' 

oo 



n=l i=l 



gn- A{q'^ 



n ■ ELi qj + 



(48) 



where we define the operator P on the right-hand side of the equation to 
order the fields with the lowest value of i to the left. From Eq. (j48|) we 
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V. - igu«t 
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Figure 3: (a) Feynman rules for eikonal lines in the amplitude and its complex 
conjugate, (b) A general contribution to a parton distribution. 



can read off the Feynman rules for the expansion of the ordered exponential 
HI] . They are illustrated in Fig. El 

The denominators n ■ qj + ie are represented by double lines, which we 
shall refer to as "eikonal" lines. These lines attach to gluon propagators via 
a vertex proportional to —ign^. Fig. El^a) shows the formal Feynman rules 
for eikonal lines and vertices. In Fig. Efb), we show a general contribution 
to fq/q, as defined by Eq. (jlUI) . 

The positions of all the explicit fields in Eq. ()46|) differ only in their plus 
components. As a result, minus and transverse momenta are integrated over. 
(They may thought of as flowing freely through the eikonal line.) The plus 
momentum flowing out of vertex 1 and into vertex 2, however, is fixed to be 
^P"*". No plus momentum flows across the cut eikonal line in the figure. Fig. 
m shows the one loop corrections to fq/q{i)- 

To be explicit. Fig. El^b) is given in n dimensions by 

— - / -^iV„^(g)27r5(g2)- \- X 

47r3y (27r)" ' {p - qf + le 

X tT{{i> - i){-igini>i'')] {igii'n^) (49) 

u ■ q — le 

where N^p is the polarization tensor of the gluon. By applying minimal sub- 
traction to Eq. and the similar forms for the other diagrams in Fig. |31 
we easily verify Eq. (j23|l for fq/q. Gluon distributions are calculated pertur- 
batively in a similar manner. We will need the concept of eikonal lines again. 
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(a) 





(d) 

Figure 4: One loop corrections to quark distribution, Eq. (pSj). 



when we discuss the proof of factorization in gauge theories. 



4.3 Renormalization 

The operator products in the definitions (pSj) and ()44|) require renormahza- 
tion, as discussed in Ref. 0. We choose to renormahze using the MS scheme. 
Of course, renormahzation introduces a dependence on the renormahzation 
scale fi. The renormalization group equation for the fa/A is the Altarelli- 
Parisi equation (jHHj) . A complete derivation of this result may be found in 
Ref. [2]. 

The one-loop result, Eq. (j2SI), can actually be understood without look- 
ing at the details of the calculation. At order as, one has simple one loop 
diagrams that contain an ultraviolet divergence that arises from the operator 
product, but also contain an infrared divergence that arises because we have 
massless, on-shell partons as incoming particles. The transverse momentum 
integral is zero, due to a cancellation of infrared and ultraviolet poles, which 
we may exhibit separately: 



j2-2e 



(2vr; 



2-2t 



1 

47r 



1 

^uv 



(50) 
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In this way, we obtain 

fa/bi^; e) = - + (— - — I - Pa/liO - counterterm + Oia^). 

(51) 

The coefficient of l/e^y is the 'anomalous dimension' that appears in the 
renormahzation group equation, that is, the Aharelh-Parisi kernel. Following 
the MS renormahzation scheme, we use the counter term to cancel l/e^y 
term. This leaves the infrared 1/e, which is not removed by renormahzation, 



(52) 



4.4 Relation to Structure Functions 

Let us now consider the relation of the parton distribution functions to the 
structure functions measured in deeply inelastic lepton scattering. If we use 
the definition of parton distribution functions given above, then the structure 
function F2 is given by the factorization equation At the Born level, 
the hard scattering function is simply zero for gluons and the quark charge 
squared, e|, times a delta function for quarks. Thus the formula for F2 takes 
the form 



X 



^F2{x, Q)= e] fj/A{x, /i) + ^? / T" ~ ( 7' ~ ) 



+ 0{al 



(53) 



The sums over j run over all flavors of quarks and antiquarks. Gluons do 
not contribute at the Born level, but they do at order as, through virtual 
quark- ant iquark pairs. The hard scattering coefficients Cjb can be obtained 
by calculating (at order as) deeply inelastic scattering from on-shell massless 
partons, then removing the infrared divergences according to the scheme 
discussed in Sect. |21 

The explicit form of the perturbative coefficients Cjb is j3] 



4 



Cjk{z, 1) — 6jk - 



11 + 



2 1 



In 



3 1 



41- z 



3 



C,,iz,l) 



2 \2 L 



In 



1- z 



+ 1 



3zil-z)\, (54) 
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where the plus subscript to the bracket in the first equation denotes a sub- 
traction that regulates the z ^ 1 singularity, 

dz [C{z)]+ h{z) = [ dz [C{z)]+ h{z) e{z > x) 

Jo 

= dzC{z) ^^h{z)eiz>x)-hil)y (55) 
4.5 Other Parton Distributions 

The definitions and ()44|1 are the most natural for many purposes. They 
are not, however, unique. Indeed, any function gb/A{y), which can be related 
to fa/A{x) by convolution with ultraviolet functions Dab{x/y, Q/ fi) in a form 
hke 

9a/A{x) = ^ (dy/y) Dab{x/y,Q/fi,as{fi)) fb/A{y) , (56) 
b 

is an acceptable parton distribution [T2] . The hard scattering functions calcu- 
lated with the distributions gt/A will differ from those calculated with /^/^(x), 
but this difference will itself be calculable from the functions Dab as a power 
series in as{Q). 

The most widely used parton distribution of this type is based on deeply 
inelastic scattering, and may be called the DIS definition. The definition is 

ff/l{x,^^) = f,/Aix,^^) + J2fJfwAi^,^^)^c,,(^^^ (57) 

for quarks or antiquarks of fiavor j. Comparing this definition with Eq. (|53|). 
we see that 

x-'F2ix,Q) = Y,e] ff/lix,Q) + Oial). (58) 
j 

That is, we adjust the definition so that the order correction to deeply 
inelastic scattering vanishes when fi = Q. It is not so clear what one should 
do with the gluon distribution in the DIS scheme. One choice ^H] is 

/,%'(^,/^) = h/A{x,f^) - f /VA(e,/i)^C7,, 1^ + 0{a^). 

(59) 
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This has the virtue that it preserves the momentum sum rule that is obeyed 
by the MS parton distributions 



If one wishes to use parton distribution functions with the DIS definition, 
then one must modify the hard scattering function for the process under 
consideration. One should combine Eqs. (j52j) and (|59|) to get the DIS distri- 
butions of a parton in a parton, then use these distributions in the derivation 
in Sect. |21 

It should be noted that there is some confusion in the literature concern- 
ing the term +1 that follows the logarithm in Cjg in Eq. (j54p . The form 
quoted is the original result of Ref. translated from moment-space to 
z-space. In the calculation with incoming gluons, one normally averages over 
polarizations of the incoming gluons instead of using a fixed polarization. 
This means that one sums over polarizations and divides by the number of 
spin states of a gluon in 4 — 2e dimensions, namely 2 — 2e. If, instead, one 
divides by 2 only, one obtains the result (jSH) without the +1, which may 
be found in Ref. ^3]. This does no harm if, as in the case of Ref. one 
wants to express the cross section for a second hard process in terms of DIS 
parton distribution functions and if one consistently divides by 2 instead of 
2 — 2e in both processes. However, it is not correct if one wants to relate 
the DIS structure functions to MS parton distribution functions, defined as 
hadron matrix elements of the appropriate operators, renormalized by MS 
subtraction. 

5 Factorization for 0^ Theory 

In this and the next section, we study the factorization theorem in a 0^ 
theory for n < 6 space-time dimensions. First we show how the factorization 
theorem comes about for one-loop corrections in deeply inelastic scattering, 
and compare the field theory to the parton model. In the next section, we 
will present a reasonably complete but compact derivation of the factorization 
theorem in deeply inelastic scattering to all orders of perturbation theory. 

The scalar theory allows us to study these issues in a simplified but highly 
nontrivial context. As emphasized above, the purpose of the factorization 
theorems is to separate long-distance behavior in perturbation theory. In the 




(60) 
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scalar theory, as we shall see, this behavior is associated with partons that 
are collinear to the observed hadrons. The organization of such "collinear 
divergences" is central to factorization in all field theories, but in gauge the- 
ories they are joined by "soft" partons, associated with infrared divergences. 
Indeed, the basic problem in gauge theories is to show how that infrared or 
"soft" divergences cancel (see Sect. EI). In 0^ theory the infrared problem is 
absent, so that studying this theory allows us to study the basic physics of 
factorization in the simplest possible setting. 
The Lagrangian is 

£ = 1 {d(f)f - im202 - -g{fih^/Any/'^(j)^ + counterterms . (61) 
2 2 6 

We will use, where necessary, dimensional regularization, with space-time 
dimension n = 6 — 2e. It is worth recalling that at n = 6 the theory is 
renormalizable, while for n < 6 it is superrenormalizable. We shall not 
concern ourselves with the theory for n > 6 where it is nonrenormalizable by 
power-counting, /i is a mass which enables us to keep g dimensionless as we 
vary n. We will renormalize the theory with the MS prescription. We use 
the factor {^'^e^ /AttY^"^ rather than the more conventional fx^, so that we can 
implement MS renormalization as pure pole counterterms. (For convenience, 
we will define the h(f) counterterm that renormalizes the tadpole graphs by 
requiring the sum of the tadpoles and their counterterm to vanish.) We 

define 

fi = /iy/eT/47r. (62) 

5.1 Deeply inelastic scattering 

Our model for deeply inelastic scattering consists of the exchange of a weakly 
interacting boson. A, not included in the Lagrangian (j6T| . This is illustrated 
diagrammatically in the same way as for QCD, in Fig. El The weak boson 
couples to the field through an interaction proportional to ^Acp"^. There is 
then a single structure function which we define by 

F{x,Q) = d'ye''^-y{p\j{y)jmp), (63) 

where j = \4>^- The momentum transfer is g'^, and the usual scalar variables 
are defined by = —q^ and x = Q'^/2p ■ q, with the momentum of the 
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Figure 5: Deeply inelastic scattering. 



target. We will investigate the structure function in the Bjorken limit of 
large Q with x fixed, and our calculations will be for the case that the state 
\p) is a single particle (with non-zero mass, as given in Eq. ()61|)). 

When Q is large, each graph for the structure function behaves like a 
polynomial of ln(Q/m) plus corrections that are nonleading by a power of 
Q. Factorization is possible because only a limited set of momentum regions 
of the space of loop and final state phase space momenta contribute to the 
leading power. First we will explain the power counting arguments that 
determine these "leading regions" , and how they are related to the physical 
arguments of the parton model. 

The tree graph for the structure function is easy to calculate. It is 



The one-loop "cut diagrams" (as defined in Sect. El above) which contribute 
to F are given in Fig. IHl 

Each of these diagrams illustrates a different bit of the physics, so we 
shall treat them in turn, starting with the "ladder" correction. Fig. IHl^a). 

5.2 Ladder Graph and its Leading Regions 

The Feynman integral for the cut diagram Fig. Efa) is 



Although for nonzero m this integral is finite, it will prove convenient to 
retain the dimensional regularization, in order to display some very important 
dimension-dependent features of the Q ^ oo limit. 

Equation ()65|) is calculated conveniently in terms of light-cone coordi- 
nates. Without loss of generality, we may choose the external momenta, 



Fo = Q^S{2p-q + q^) = S{x-l). 



(64) 




(65) 
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Figure 6: One-loop corrections to deeply inelastic scattering. For graphs (b), 
(d) and (e), we also have the hermitian conjugate graphs. 



and to be = (p"*", m^/2p+, 0_l) and = (— Q^/2xp+, 0_l). Notice 
that this formula for q'^ corresponds to a slight change in the definition of x, 
which we now define by Q^/2p ■ q = x/{l — xm? /Q"^). At leading power in 
there is no difference, but at finite energy our formulas will be simplified 
by this choice. 

The 5-functions in fl65p can be used to perform the k± and k' integrals. 
Then if we set ^ = k'^ jp^ ^ we find 



2(a) 



647r^ \e"i ix^x{\ — x"^ 



r(2 



2 \ 1-e 

m X 



X 



1-e 



x[(^max-0(^-^min)] 



(66) 



where the limits S,rain and ^^ax are given by 



1 + X ^1 — X 



Am'^x 



[1 — x){Q'^ + m?x) 



(67) 



In this form, we can look for the leading regions of the ladder corrections. To 
do this, it is simplest to set the mass to zero, find the leading regions, and 
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then check back as to whether we must reincorporate the mass in the actual 
calculation. So, to lowest order in A = rn? /Q"^, becomes 



To interpret this expression, we must distinguish between the renormalizable 
= 6, e = 0) and superrenormalizable (ra < 6, e > 0) cases. 
In the super-renormalizable case, (e > 0), the leading-power contribution 
[Q/n)^ comes from near the endpoint ^ = x{l + A). The bulk of the integra- 
tion region, where ^ — x = 0(1) is suppressed by a power of Q. The integral 
is power divergent when m = 0, and clearly we cannot neglect the mass. 

Now consider the renormalizable case, n = 6. When we set e = 0, Eq. ()(i8j] 
has leading power (Q°) contributions from both the region ^ — x near zero, 
where, as above, the mass may not be neglected, and the region ^ — x = 0(1), 
where it may. In the former region, the integral is logarithmically divergent 
for zero mass, but since the nonzero mass acts as a cutoff, the two regions 
^ ~ x and ^ — X = 0(1) should be thought of as giving contributions of 
essentially equal importance. We now interpret these dimension-dependent 
leading regions. 

5.3 Collinear and Ultraviolet Leading Regions; the 
Parton Model 

To see the physical content of the leading regions identified above, it is useful 
to relate the variable ^ in to the momentum /c^, by the relations 





l-e 



(68) 





(70) 
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Figure 7: Space-time structure of collinear interaction. 



We emphasize that this expression is accurate to leading power in the region 
IQ^ = 0(A), which is sufficient to give the full leading power for n < 6, 
although not for n = 6, where larger k±_ also contribute. 

Now let us choose a frame in which p'^ is of order Q. When ^ — x, the 
components of are of order (Q, — x)Q, Qy/^ — x), and at its lower hmit, 
^ — X is of order m? /Q"^. Hence, in the region that gives the sensitivity to m, 
fcx is small, and is ultrarelativistic and represents a particle moving nearly 
collinear to the incoming momentum, p^. In addition, the on-shell line, of 
momentum p^ — k^, is nearly collinear to the incoming line as well. In fact, 
when m and k± are both zero, k^ is also on the mass shell. The energy 
deficit necessary to put both the momenta k^ and p^ — k^ on shell is of order 
k±^ /Q in this frame. Thus, in this frame, the intermediate state represented 
by the Feynman diagram lives a time of order Q/fc^^, which diverges in the 
collinear limit. The space-time picture for such a process is illustrated in 
Fig. [71 and we see a close relation to the parton model, as discussed in Sect. 
11.41 which depends on the time dilation of partonic states. Partonic states 
whose energy deficit is much greater than m in the chosen frame correspond 
to ^ — X of order unity, and do not contribute at leading twist. Thus here, 
as in the parton model, there is a clear separation between long-lived, time- 
dilated states which contribute to the distribution of partons from which the 
scattering occurs, and the hard scattering itself, which occurs on a short time 
scale. 

From this discussion, the collinear region, which is the only leading region 
when n is less than 6, is naturally described in parton model language. 
When n = 6, the collinear region remains leading. In addition, however, 
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all scales between kj_ = m and kj_ = Q contribute at leading power, and there 
is no natural gap between long- and short- distance interactions. When ^ — x 
is order unity, fc^ is separated from p'^ by a finite angle, and corresponds 
to a short-lived intermediate state, where [p — k)"^ ^ m^. This leading 
region, which is best described as "ultraviolet" , is not naturally described by 
the parton model. But, in an asymptotically free theory (as (0^)g is), such 
short-lived states may still be treated perturbatively. We shall see how to do 
this below. 

In summary, the ladder diagram shows two important features: a strong 
correspondence with the parton model from leading coUinear regions for 
both superrenormalizable and renormalizable theories and, for the renor- 
malizable theory only, leading ultraviolet contributions, not present in the 
parton model. 



5.4 Parton distribution functions and parton model 

We shall now freely generalize the results for the one loop ladder diagram. 
Indeed, as we shall see in Sect. 13 some of the dominant contributions to 
the structure function arise from (two-particle-reducible) graphs of the form 
of Fig. |H| A single parton of momentum k'^ comes out of the hadron and 
undergoes a collision in the Born approximation. If we temporarily neglect 
all other contributions, we find that 

F{x,Q) = I ^^mp)H{k,q) + 0{l/Qn, (71) 

where $ represents the hadronic factor in the diagram and H the hard scat- 
tering (multiplied by the factor of Q'^/2tt in the definition of the structure 
function): 

H{k,q) = Q^6iik + Qy-m^). (72) 

For 0^ with ri < 6, as in the parton model, the parton momentum k'^ 
is nearly collinear to the hadron momentum p^. This implies that we can 
neglect m and the minus and transverse components of k'^ in the hard scat- 
tering, so that we can write 

H = H{x/^,Q) = S{ax-l), (73) 

and hence 



F{x,Q)= / d^ 







,dk d^fc^ 



6{^/x-l) + 0{l/Q^). (74) 
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Figure 8: Dominant graphs for deeply inelastic scattering in parton model. 

Here we define ^ = jp^ . The limits on the ^ integral are to 1, since 
the final state must have positive energy. We therefore define the parton 
distribution function (or number density): 



d/t-d^fci 



(2vr)^ 



$(fc,p) 



(2^ 



$(A;,p)5(^j9+/A;^ 



(75) 



With this definition ()74|) becomes 

F(x,Q) = J^f(^)H{x/^,Q) + 0{l/Q'^) 



f{x) + Oil/Q^) (n<6). 



(76) 



As we shall see in the next section, the factorization theorem is also true in 
the renormalizable theory, 



F{x,Q) 



fiOHix/^,Q) + Oil/Q^) 



[n 



6) 



(77) 



where now H is nontrivial. The dominant processes that contribute are 
illustrated by Fig. El which generalizes the parton model only to the extent 
of having more than just the Born graph for the hard scattering. These 
processes first involve interactions within the hadron that take place over 
a long time scale before the interaction with the virtual photon. Then one 
parton out of the hadron interacts over a relatively short time scale. 
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Figure 9: Dominant regions for deeply inelastic scattering in (<^^)g theory. 



\ 0V\ A 



Ct) tt>' tc) 

Figure 10: Low-order graphs for parton distribution in 0'^ theory. 

We now note that Eq. (f75|) can be expressed in operator form as 



/(O 



2tx 



-i^p+y 



{p\4>{0,y-,0^)m \p). 



(78) 



This is the definition which we use for all n < 6. Of course in the renormal- 
izable theory n = 6 renormalization will be necessary [H]. The definition ()78|) 
is precisely the analog for 0^ theory of those we gave in Sect. |3]for QCD. It 
involves an integral over a bilocal operator along a light-like direction. The 
graphs for /(,^) up to one- loop order are shown in Fig. Feynman rules 
are the same as for the gauge theory, but without eikonal lines. 

It is natural to interpret f{x)dx as the number of partons with fractional 
momenta between x and x + dx. This interpretation is justified by the use 
of light front quantization .10., as we saw in Sect. |3] Note that although 
the definition picks out a particular direction as special to the problem, it is 
invariant under boosts parallel to this direction. 
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The ladder graph, Fig. ITUT c) gives 



(79) 



The 5-functions may be used to perform the k'^ and k integrals, after which 
we obtain 

^^-647r3 r(2-e) ], [fc/ + ^^(1 - ^ + ^^^^ ^ ^ 



which matches Eq. ()70|) in the fc^ limit. That is, we have constructed 
the parton distribution to look like the structure function at low transverse 
momentum. The significance of this fact will become clear below. 

For n < 6, Eq. (jHDI) is the same as the full leading structure function (fTUjl . 
and it exemplifies the validity of the parton model in a super-renormalizable 
theory. When n = 6, however, there is a logarithmic ultraviolet divergence 
from large fc^ in (jHUj) . So, in the renormalizable theory we must renormalize 
/(^). (Since / is a theoretical construct defined to make treatments of high- 
energy behavior simple and convenient, we are entitled to change its definition 
if that is useful; in particular, we are allowed to include renormalization in 
its definition.) If we use the MS scheme, then the renormalized value of fc 
for nonzero mass is: 



while for zero mass it is (compare Eq. 



(81) 



i?[/e] = -^e(l-Oj. (82) 

Now let us see what this means in the calculation of the hard part, as in Sect. 
121 To calculate the hard part, we expand Eq. (ff7|l in powers of g^, as in Eq. 
(121, and solve for H^^\x/i,Q). There is some question about what to do 
with the higher-twist terms, proportional to powers of m/Q. The simplest 
method is to simply define 

H^'^ (x/e, Q) = (x/e, g) - /(^)] _o , (83) 

Comparison of Eqs. ()70|1 and (j80|l shows that the low k± region, which is 
the only leading region which is sensitive to the mass, cancels between F^^^ 
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and f^^\ at the level of integrands. Thus, for the combination on the right 
hand side of Eq. (|HH|). it is permissible to set the mass to zero. It is thus 
practical to set the mass to zero at the very beginning. It should be kept in 
mind, however, that this is a matter of calculational convenience, rather than 
principle. The factorization theorem allows us to calculate mass-insensitive 
quantities whatever the masses we choose, since all sensitivity to these masses 
will be factored into the parton distributions. 

Now let us return to the remaining diagrams in Fig. ^ treating first the 
"final state" interactions. Fig. IHl^b) and (c). 



5.5 Final state interactions 

The graphs of Fig. Efb) and (c) have a self-energy correction on the outgoing 
line, the final state cut either passing through the self energy or not. As we 
will show, these graphs have contributions that are sensitive to low virtual- 
ities and long distances. However, they are not of the parton model form, 
and do not naturally group themselves into the parton distribution for the 
incoming hadron. We will see, however, that there is a cancellation between 
the two graphs such that they are either higher twist (n < 6), or may be 
absorbed into the one- loop hard part (n = 6). 

The self energy graphs give simply the lowest order graph, 6{x — 1), times 
the one- loop contribution to the residue of the propagator pole: 



2(b) 



5{x - 1) 



1287r 



^0 



d^ / d/c_L A;_L 



-f- counterterm 



We may derive this expression in either of two ways. One way is to combine 
the denominators of the two propagators in the loop by a Feynman parameter 
before performing the and k~ integrals. Then z is the Feynman parameter. 
Alternatively, we may first use contour integration to perform the k~ integral. 
Then we get (j84j) by writing = z {p^ + The integral is the same by 
either derivation. But the second method shows that we may interpret z as 
a fractional momentum carried by one of the internal lines. Since we will be 
concerned with the low k±'^ region, while the counterterm, if computed with 
MS renormalization, is governed by the k± ^ oo behavior of the integrand, 
we do not write the counterterm explicitly. 
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There is clearly a significant contribution in (j84j) from small k^, where 
the mass m is not negligible. The cut self-energy graph, in Fig. El^c), will also 
contribute in this region. Now the region of low k± represents the effect of 
interactions that happen long after the scattering off the virtual photon, and 
it is reasonable to expect that interactions happening at late times cancel, 
since the scattering off the virtual photon involves a large momentum transfer 
Q and therefore should take place over a short time-scale. However, the uncut 
self-energy graph only contributes when x is exactly equal to 1, while the cut 
self energy graph has no (5-function and thus contributes at all values of x. 

This mismatch is resolved when we recognize that we should treat the 
values of the graphs as distributions rather than as ordinary functions of x. 
That is, we consider them always to be integrated with a smooth test func- 
tion. Mathematically, this is necessary to define the J-function. Physically, 
the test function corresponds to an averaging with the resolution of the ap- 
paratus that measures the momentum of the lepton that is implicitly at the 
other end of the virtual photon. After this averaging, a measurement of the 
lepton momentum does not distinguish the situation where a single quark 
goes into the final state from the situation where the quark splits into two. 

We therefore consider an average of the structure function F{x) with a 
smooth function t{x): 



{t,F) = J dxt{x)F{x). (85) 

Then the contribution of the self-energy graph is 

(t,F2(fe)) = t(l)F2(,) . (86) 
Next we compute the cut graph. Fig. Efc). Its value is 



(t,F2(c)) = l^/d^^/ dk^'k^' I j dk- 



t{x) S{k'^ — m^) 5{{p + q — ky — m^) 
[{p + ~ ''^^]^ 

To make this correspond with the form of ()84|1 . we define z = k'^ /{p'^ + q^), 
and then use the 5-functions to do the k~ and x integrals. After some algebra, 
and after the neglect of terms suppressed by a power of Q, we find 

, r. ^ 9^ /"\ P 2 , 2 Z)t(x)x'^ 



40 



where the Bjorken variable x satisfies 



X 



1 + 



Q^z{l - z) 
We now add the two diagrams to obtain: 

{t-, F2{b) + -F2(c)) 



9' 



dz 



1287r3 J, J, 
+ counterterm. 



dk^^k^^[t{x)x^ -t{l)] 



(89) 



z{l-z) 



[k^' + m^{l- z + z^)Y 



In the region k± <^ Q, x is close to one, and there is 
integrand of Eq. ()9()j) . The cancellation fails when z is 
but the contribution of that region is suppressed by a 
therefore permitted to set m = in the calculations 
which a calculation (with dimensional regularization to 
divergences that now appear in each individual graph 
easier. 



(90) 

a cancellation in the 
close to zero or one, 
power of Q. We are 
of the graphs, after 
regulate the infrared 
at k± = 0) is much 



5.6 Vertex correction 



Finally, we consider the vertex correction Fig. IHfd). It has the value 



6{x - 1) 



6-2e 



(2vr; 



-2e; 



[m^ — A;^] [m^ — {p + ky] [m? — [p + k + g)^] 
+ counterterm 

,2 



-6{x-l) 
X In 



9 



dai 



1— ai 



dag 



647r3 

m^(l - ai - q;2 - ("i + "2)^) + <5^a2tti 
. J' 



(91) 



where we work in = 6 — 2e space-time dimensions to regulate the ultraviolet 
divergence. When Q 00, we can clearly neglect the mass, so that we have 
(at e = 0) 



-5{x - I) 



9' 



1287r3 



o 



1 



(92) 
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F2[d) is higher twist for e > 0. 

The graph Fig. IHl^e) is related to Fig. IHl^d) by moving the final state cut 
so that it cuts the inner lines of the loop. We will not calculate it explicitly. 
But when that is done, the quark mass can be neglected, just as for the uncut 
vertex. 

In summary, the only diagram from Fig. El which corresponds to the par- 
ton distribution is the ladder diagram. Fig. E^a). Non- ladder diagrams are 
either higher twist, or contribute only to the hard part (renormalizable case). 
These results are consistent with the structure of Fig. |H1 and Fig. El which 
show the structure of regions which give leading regions for n < 6 and n = 6, 
respectively. As we shall show in the next section, it is this structure which 
enables us to prove that the parton distributions Eq. (f7S|) absorb the com- 
plete long-distance dependence of the structure function. 

6 Subtraction Method 

To establish a factorization theorem one must first find the leading regions 
for a general graph. We will see how to do this in Sect. [71 The result, for 
deeply inelastic scattering in a nongauge theory, has been summarized by 
the graphical picture in Fig. El and it corresponds closely to our detailed 
examination of the order g"^ graphs. It can be converted to a factorization 
formula if one takes sufficient care to see that overlaps between different 
leading regions of momentum space do not matter. 

An approach that makes this process clear is due to Zimmermann [2|. 
To treat the operator product expansion (OPE), he generalized the methods 
of Bogoliubov, Hepp, Parasiuk, and Zimmermann (BPHZ) |2l El that were 
used to renormalize Feynman graphs. (Although the original formulation 
was for completely massive theories with zero momentum subtractions, it 
can be generalized to use dimensional continuation with minimal subtraction 
[inj. This allows gauge theories to be treated simply.) In the case of deeply 
inelastic scattering a very transparent reformulation can be made in a kind 
of Bethe-Salpeter formalism [17j, although it is not clear that in the case of 
a gauge theory the treatments in the literature are complete. In this section, 
we will explain these ideas in their simplest form. 

There are two parts to a complete discussion: the first to obtain the 
factorization, and the second to interface this with the renormalization. We 
will treat only the first part completely. In (0^)g theory, renormalization 
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is a relatively trivial affair. Moreover, if we regulate dimensionally, with e 
just slightly positive, one can choose to treat as the leading terms not only 
contributions that are of order (times logarithms) as Q — oo, but also 
those terms that are of order Q to a negative power that is of order e. The 
remainder terms are down by a full power of Q^, and can be identified as 
"higher twist" . In this way one has the same structure for the factorization, 
without the added complications of renormalization. 

Zimmermann's approach is to subtract out from graphs their leading 
behavior as Q — oo. This is a simple generalization of the renormalization 
procedure that subtracts out the divergences of graphs. From the structure 
function Q) one thereby obtains the remainder -FRem(a;, Q), which forms 
the higher twist contributions. The leading twist terms are F — F^em- It is a 
simple algebraic proof to show that F — -Faem has the factorized form f * H, 
with / being the parton distribution we have defined earlier, and with 
denoting the convolution in Eq. (f77|). 

6.1 Bethe-Salpeter decomposition 

In the graphical depiction of a leading region. Fig. El exactly one line on 
each side of the final state cut connects the collinear part and the ultraviolet 
part. So it is useful to decompose amplitudes into two-particle- irreducible 
components. This will lead to a Bethe-Salpeter formalism. Consider, for 
example, the two-rung ladder graph. Fig. ^2 for deeply inelastic scattering 
off a composite particle. We can symbolize it as 

Fig. [n]= 7s X 7 X 7 X (93) 

Here 7^ represents the graph that is two-particle-irreducible in the vertical 
channel and is attached to the initial state particle, 7 represents a rung, and 
7/i represents the two-particle-irreducible graph where the virtual photon 
attaches. It is necessary to specify where the propagators on the sides of 
the ladder belong. We include them in the component just below. Thus 7^ 
and 7 have two propagators on their upper external lines. The purpose of 
having a composite particle for the initial state is to give an example with a 
non-trivial 7s, as in QCD with a hadronic initial state. The vertex joining 
the initial particle is a bound-state wave function. 
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Figure 11: Example of ladder graph with several rungs. 
We now decompose the complete structure function as 

oo 

F = Gs Gh 

N=0 

= G^Y^^G,. (94) 

Here Gg is the sum of all two-particle- irreducible graphs attached to the initial 
state particle, G^ is the sum of all two-particle-irreducible graphs coupling to 
the virtual photon, and Gr is the sum of all graphs for a rung of the ladder. 
Thus Gr is the sum of all two-particle-irreducible graphs with two upper lines 
and two lower lines, multiplied by full propagators for the upper lines. 

The second line of Eq. fl94|) has the inverse of 1 — G^, and it clearly 
suggests a kind of operator or matrix formalism. Indeed, if we make explicit 
the external momenta of two ladder graphs, 7i(A;,/) and 72(A;,/), then their 
product is 

(7i72)(fc,0 = J ^^7i(^,A:')72(A:',/). (95) 

The rung graphs can thus be treated as matrices whose indices have a con- 
tinuous instead of a discrete range of values, while Gg and Gh can be treated 
as row- and column-vectors. 

In the case that the initial hadron is a single parton, as in the low order 
examples in Sect. El the soft part Gg is trivial: Gg = 1, where '1' represents 
the unit matrix. 
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Figure 12: Hard scattering coefficient from Fig. IHfa). 

6.2 Extraction of higher twist remainder 

We can now symbolize the operations used to extract the contribution of a 
graph to the hard scattering coefficient. Consider the example that lead to 
Eq. (j83|) . We took the original graph and subtracted the contribution of the 
graph to f^^^H^^\ where H^^'^ is the lowest order hard part. Then we took 
the large Q asymptote of the result, by setting all the masses to zero. 

We represent this in a graphical form in Fig. ^1 There, the wavy line 
represents the operation of short circuiting the minus and transverse compo- 
nents of the loop momentum coming up from below, and of setting all masses 
above the line to zero. Symbolically, we write this as: 

Contribution of Fig. Efa) to H'^^^ = P'j^h — PlPlh 

= P7(l-P)7,, (96) 

where the operator P is defined by 

P(fc, /) = (27r)6- 5{k^ - /+) 5{n 6'-^{h) ^^^^ 
X (Set masses to zero in the part of the graph above P). 

In Eq. fl96|) we have ignored the need for renormalization that occurs if e = 0. 
Either we can assume that we are only making the argument when e is slightly 
positive, or assume that all necessary renormalization is implicitly performed 
by minimal subtraction. 

Fig. infa) gives two contributions to the factorization: a contribution to 
the one-loop hard part H^^^ given in Eq. ()96|) or (jHHj) . and a contribution to 
y (1)^(0) ^j^g second of these we picture in Fig. [121 and symbolize as 

iPlh. (98) 
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Figure 13: Contribution of Fig. a) to f^^H^^\ 

Thus we can write the remainder for Fig. IHfa), after subtracting its leading 
twist contribution, as 

Rem(Fig.inta)) = 77?,, - - P7(l - P)7fe 

= (1-P)7(1-P)7;,. (99) 

Clearly the operator 1 — P subtracts out the leading behavior. 

In general, we can write the remainder for the complete structure function 

as 

00 

FKem = J]a[(i-p)a]^(i-p)G^ 

N=0 

This formula is valid without renormalization, even at e = 0. In the first 
place, renormalization of the interactions can be done inside the 7's. This 
is because there is nesting but no overlap between, on the one hand, the 
graphs to which the operation 1 — P is applied and, on the other hand, 
the vertex and self-energy graphs for which counterterms are needed in the 
Lagrangian of the theory. Further divergences occur because of the extraction 
of the asymptotic behavior, and these give rise to the need to renormalize the 
parton distribution. But the regions that give rise to such divergences are 
of the form where lines in some lower part of a graph are coUinear relative 
to lines in the upper part. All such regions are canceled in Eq. ()100j) since 
to the operator 1 — P they behave just like the regions that give the leading 
twist behavior of the structure function. 
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6.3 Factorization 



It is now almost trivial to prove factorization for the leading twist part of 
the structure function, which is 

F — -FRem = Gs ^ _ Q — Gs{l — P) ^_ Q ^-[^ _ p'j (101) 

Simple manipulations give 

= f*H. (102) 
We now have an explicit formula for the hard scattering coefficient: 

while the parton distribution / satisfies 

/xP = G,^^P. (104) 

One somewhat unconventional feature of our procedure is that not only 
do we define P to set to zero the minus and transverse components of the 
momenta going into the subgraph above it, but we also define it to set masses 
to zero. Setting the minus and transverse momenta to zero while preserving 
the plus component is exactly the appropriate generalization of BPH(Z) zero- 
momentum subtractions to the present situation. Setting the masses to zero 
as well is a convenient way of extracting the asymptotic large-Q behavior 
of a graph, as we saw in our explicit calculations. Moreover, particularly in 
QCD, it greatly simplifies calculations if one works with a purely zero-mass 
theory. Of course, setting masses to zero gives infrared divergences in all but 
purely ultraviolet quantities. The momentum-space regions that give the di- 
vergences associated with the structure function all have the same form as 
the leading regions for large Q, Fig. El so that the 1 — P factors in Eq. p04p 
kill all these divergences. Note that, just as with Zimmermann's methods, 
the P operator can be applied at the level of integrands. In practical calcula- 
tions, dimensional continuation serves as both an infrared and an ultraviolet 
regulator. 
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In the one-loop example of Sect. El the external hadron is a parton, so 
that in Gg = S{x — 1) in (|l(J4p . At one loop, Gr corresponds exactly to 
fc, Eq- (jHUI) . This expression, and the distribution / as a whole in p04p 
is still unrenormalized, and contains ultraviolet divergences. These may be 
removed by minimal subtraction, as in Eq. (jHT|) at one loop, or as discussed 
more generally in Ref. 0. We should mention, however, that it would be 
advantageous to have a subtraction procedure which combined factorization 
and renormalization into a single operation. The particular procedure out- 
lined by Zimmermann [2] does this, but is not immediately applicable when 
all particles are massless. Duncan and Furmanski jTH] have discussed some 
of these issues at length. 

6.4 Factorization for Inclusive Annihilation in (0^)g 

It is easy to generalize the general arguments of this and the previous section 
to other processes, such as those listed in the introduction. An important 
example, is the cross section in 0^ theory that is analogous to one-particle 
inclusive annihilation in e~^e~ annihilation, that was discussed in Sect. 11.21 
In the scalar theory, the structure function for this process is 

Dix,Q) = ^ [ d%e''^-yY,{0\j{y)\HX){HX\j{0)\0), (105) 

which is exactly analogous to the QCD version, Eq. ©. 

It is relatively easy to check that the leading regions for this process have 
a form that generalizes Fig. 121 for deeply inelastic scattering, that is, they 
have the form of Fig. This was shown in Ref. (for the case of a 
non-gauge theory). 

An example is given by the ladder graph of Fig. ^1 We must integrate 
over all values of the momentum {k — p)^- When {k — p)'^ is collinear to 
pf^, the line k has low virtuality. Then in the overall center-of-mass, the 
remaining particle q — k has large energy, approximately Q/2, and is moving 
in the opposite direction to the first two particles. We therefore consider the 
lines p, k — p and k as forming the jet Ja in Fig. and q — k together with 
the vertex where the 'virtual photon' attaches as forming the hard part H. 
When [k — p)'^ has transverse momentum of order Q, we put both k and 
k ~ p into the hard part. 

In a non-gauge theory, these two regions are the only significant ones, 
together with a region that interpolates between them. As we shall see in 
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Sect. 13 this statement generalizes to all orders of perturbation theory. In 
a gauge theory, like QCD, all kinds of complication arise because there are 
also 'leading twist' regions involving soft gluons. 



6.5 Factorization, fragmentation function 

Simple generalizations of the arguments for deeply inelastic scattering give 
the scalar factorization theorem: 

D{z, Q) = j'j H{z/C, Q) d{C) + 0(l/g"). (106) 

analogous to Eq. ((Tj). Here the fragmentation function is defined in exact 
analogy to the parton distribution. We choose axes so that the momentum 
p^^ of the detected particle is in the positive 2;-direction. Then we define: 



d{z) = ^ [dy-e'''^y-/^Y.(^\<P{0,y~,OT)\HX){HX\m\0) 

-^Dik,p). (107) 



X 



z J (27r) 

This is interpreted as the number density of hadrons if in a parton. The for- 
mulae are exactly analogous to Eqs. (ffK|l and (|7H|l for the parton distribution. 
Renormalization is needed here also. 



7 Leading Regions 

As we saw in Sects. ElandlHl the first step in constructing a complete proof of 
a factorization theorem is to derive the leading regions of momentum space 
for a graph of arbitrary order. This section begins with a brief description 
of a general approach to the long- and short-distance behavior of Feynman 
diagrams that results in a derivation of the leading regions. We apply this 
method to describe the origin of high-energy logarithms in scalar theories, 
and go on to discuss the cancellation of final state interactions, and the 
infrared finiteness of jet cross sections. 
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7.1 Mass dependence and leading regions 

Consider, then, an arbitrary Feynman integral G{p'^ / fi,m/ fi), corresponding 
to a graph G, which is a function of external momenta pj', mass m (possibly 
zero), and renormalization scale /x. Without loss of generality, we may take 
G to be dimensionless. We also assume that the invariants formed from 
different p'^ are all large, while the themselves have invariant mass of 
order m. Thus: 

Pi ■ Pj = Q'^Vij, p1 = (108) 

where Q is a high-energy scale, Q ^ m, and the rjij and Q are numbers 
of order unity. In the following, it will not be necessary to consider the rjij 
and (i dependence, and we will write G as G(Q^/yU^,m^//i^). We will be 
interested in the leading term in an expansion in powers of l/Q^. (Always 
we will allow the possibility of a polynomial in In Q multiplying the power of 
Q, in each order of perturbation theory.) 

Suppose G is the result of L loop momentum integrations acting on a 
product of / Feynman propagators, times a function N, which is a polynomial 
in the internal and external momenta. For simplicity, we absorb into the 
numerator factors associated with the internal propagators, as well as overall 
kinematic factors, etc. G may then be represented schematically as 

G{Qy^,^ mV/i') = n / d"£, N{{k,}, {p,}) n \^.y (109) 

i=i J j=i y'^j ^ ^ 

The line momenta {^^ }, of course, are functions of the £f and the . Any 
region in ^^ space which contributes to G at leading power in will be 
called a "leading region". In addition, by a "short-distance" contribution 
to p09p we will mean that we have a region of loop momenta in which 
some subset of the line momenta, {^^ }, are off-shell by at least (9(Q^); the 
short- distance contribution is the factor in (|1U9|) given by these far off-shell 
lines. Short- distance contributions are independent of masses to the leading 
power in Q^, since the integrand can usefully be expanded in powers of 
when propagators are far off-shell. A general leading region has both short- 
and long-distance contributions, the latter associated with lines which are 
nearer the mass shell. Roughly speaking, factorization is the statement that 
the cross section is a product of parton distributions, in which all the long- 
distance contributions are found, and a hard-scattering coefficient, which 
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has purely short- distance contributions. To study factorization, we must 
characterize all "long- distance" contributions. 

Our analysis depends on two observations. The first concerns the close 
relation between the high-energy and zero-mass limits. That is, if the renor- 
malization scale /i is chosen to be of 0{Q), then the two limits are equivalent 
in the function G{l,m? /Q"^). Short-distance contributions to the Q ^ oo 
limit are those involving lines for which fc| is of order . Long-distance 
contributions are parts of the £f integrations for which /c| is much less than 
Q^. If we scale all momenta down by a factor proportional to Q, then we are 
considering the m — >^ limit instead. The short- distance contributions now 
have fixed and the long-distance contributions have Feynman denomina- 
tors /c| + ie in Eq. (jl(J9|) that vanish in the m — limit. 

Note that if G is such that it only has short distance contributions, then 
the Q — > cxD limit is 0), i.e., we can just set m = 0. The QCD coupling, 
as{Q)i is an implicit argument for G, and we have already chosen to set the 
renormalization scale /i equal to Q. Thus in this case the detailed large Q 
behavior is renormalizat ion-group controlled in a simple way. 

When there are long-distance contributions to G, an expansion in powers 
of m will often fail. So to find the long-distance contributions to G, one must 
look for singularities in the m — limit. There are apparent exceptions to 
this rule, exemplified by the integral 



However, if we factor out the numerator factor m^, we are left with an integral 
that is singular like Xjvr?. This singularity is governed by the denominator. 
So what we are looking for is singularities in the m dependence in the integral 
over the denominators of G. 

Our second observation is that the integrals in ()109|) are defined in com- 
plex £f -space. As a result, it is not enough for a set of denominators to vanish 
in the integrand of ()109|) for the integral to produce a singularity at m = in 
G. We must have, in addition, a pinch of one or more of the £f integrals at 
the position of the singularity, between coalescing poles. This fact enables us 
to apply the simple but powerful analysis due originally to Landau 1^ 
on the relation of singularities in Feynman integranc^s to the singularities of 
Feynman integra/s. In the next subsection we explain the application of this 
argument. 




(110) 
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7.2 Pinch surfaces 



We begin by using Feynman parameterization to combine the denominators 
of Eq. (Unni) by 



(/-i)!n rda,/i-x:«,)n / 



N{{k,},{p,]) 



'1111 



where we have exhibited the loop-momentum dependence of the hne mo- 
menta. There is now a single denominator D{ii,aj), which is quadratic in 
loop momenta and linear in Feynman parameters. Suppose D{ii,aj) van- 
ishes for some value of loop momenta and Feynman parameters. We will 
now derive necessary conditions for this zero to produce a singularity in G. 
Then we will apply these conditions to the case m = 0. 

A pole from D = will not give a singularity in G if D can be changed 
from zero by a deformation that does not cross a pole in any one of the 
momentum or parameter contours. Consider first the parameter integrals. 
Because D is linear in the {aj}, a deformation of the aj integral will change 
D away from zero, unless 



m or a,- 







(112) 

for each line. In the first case, D is independent of aj, while in the second 
we note that aj = is an endpoint of the aj integral, away from which it 
cannot be deformed. 

Now suppose pi2j) is satisfied, and consider the momentum integrals. D 
will be independent of those loop momenta which flow only through lines 
whose Feynman parameters are zero. The contours of the remaining loop 
momenta must be pinched between singularities associated with the vanishing 
of D. Since D is a quadratic function of the remaining momenta, each 
momentum component sees only two poles in its complex plane due to the 
vanishing of D. The condition for a pinch is thus the same as the condition 
that the two zeros of the quadratic form be equal. That is, in addition to 
D = we must have dD/di^ = for all £f which flow through one or more 
on-shell lines. For each such loop momentum, the extra condition is 



0, 



;ii3) 
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where the sum goes over all lines through which the loop momentum flows. 
(Note that any line which is not on shell has Uj = 0, by Eq. (jll2p . so the 
condition ()113|) can be applied to every loop.) Together, ()112|) and ()113p are 
known as the "Landau equations" . We shall refer to any surface in momen- 
tum space on which the Landau equations are satisfied as a "pinch surface" 
of the diagram G. With each pinch surface we associate a "reduced dia- 
gram" , in which all off-shell lines are shrunk to points. By construction, the 
reduced diagram contains only those loop momenta of the original diagram 
which satisfy (jll3p with nonzero a's. 

7.3 Physical propagation 

The Landau equations are surprisingly restrictive, especially in the zero-mass 
limit. To see why, let us rederive the observation of Coleman and Norton 
1221 12H] that Eqs. ()112|) and ()113|) have an appeahng physical interpretation. 
Consider a given pinch surface. We rewrite (jll3|) on this surface as 

^(a, u,)v^ = 0, (114) 
j 

where Vj and ooj are the four-velocity and energy associated with the mo- 
mentum kj. The units of the Feynman parameters are arbitrary, so suppose 
we may, if we wish, interpret aj as the frame-independent ratio of a time to 
the energy of line j. Then each of the components of the vector {ajUjj)Vj 
has the units of a distance in space-time. It is the distance traversed in time 
(djUJj) by a free particle moving classically with velocity Vj. 

Now suppose we associate a definite position to one of the vertices 
in the reduced diagram associated with the pinch surface. Then, if line j 
attaches to the vertex at x^, x^+{aj ujj)Vj may be interpreted as the position 
of the vertex at the other end of line j. Continuing in this manner, we can 
associate with the reduced diagram a position in space-time for every one 
of its vertices, and a physical process in which free particles move between 
these points. Equation ()113|) ensures that this program can be carried out 
consistently, by requiring that in going around any closed loop we come back 
to the same position. We can use this construction as a necessary condition 
for a pinch surface. 

Finally, note that, because Eqs. ()112|) and ()113p are homogeneous in 
the a's, a rescaling of the a's leaves the Landau equations satisfied. Hence 
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the vertices in the physical picture are an indefinite distance apart, and, in 
particular, this distance may be arbitrarily large. 



7.4 Collinear and infrared pinches; power counting 

For a general diagram with arbitrary masses and external momenta the cri- 
terion of physical propagation allows a very rich analytic structure. In the 
massless limit, however, this structure actually simplifies since multiparticle 
thresholds become degenerate. The physical processes of which an isolated 
massless (but not massive) particle is kinematically capable are as follows. 
First, a massless particle of momentum p^^ may split into two (or more) 
massless particles of momenta ap^ and (1 — a)p'^, and vice- versa. This is 
the source of collinear divergences. Second, a particle may emit or absorb 
one (or more) zero energy particles. This is the source of infrared diver- 
gences. We easily check that arbitrary loops involving only collinear and 
zero-momentum particles can satisfy the Landau equations. Generally, we 
will describe a subdiagram consisting of mutually collinear particles as a "jet" 
subdiagram. Lines that have zero momentum in the massless limit we will 
call "infrared" . A jet subdiagram describes the evolution of a set of collinear 
lines, as they absorb and emit other collinear and infrared lines. 

As an example, let us consider the vertex correction Fig. IT^ . in dimen- 
sionally regulated 0^ theory. 



Vip,p') 



d"fc 1 



(27r)" [{p' — ky — VP? + ie] [(p + A;)^ — + it\ — vn? + ie) ' 

(115) 

where we assume a production process, = {p' + pY > 0. This example is 
used for illustrative purposes only. The term "leading" will refer here only 
to this diagram, and not to the behavior of the Born diagram. In fact. Fig. 
ITM a) is nonleading compared to the Born process for all n < 6. 
The Landau equation for Fig. IT^ a) is 

ai{p' - k)'' + a2{p + k)^ + aak^ = 0. (116) 

For on-shell (p^ = p'^ = m?) scattering with m 7^ 0, Eq. ()116|) has no 
solutions at all. Note that this is the case even though there is a singular 
surface illustrated in Fig. IW b). with 

(p + kf = {p' - kf = m^, k^ < 0. (117) 
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Figure 16: (a) Vertex correction, (b) Reduced diagram corresponding to Eq. 
((TTTj), (c) Infrared reduced diagram, (d) coUinear reduced diagram. 

This singular surface corresponds to the production of two particles, followed 
by a subsequent spacelike scattering. Although such a process is kinemati- 
cally possible, it clearly cannot correspond to physical propagation, because 
the two particles produced at vertex 1 propagate in different directions, and 
would therefore not be able to meet at vertex 2 to scatter again. 

Now let us consider the case that m = 0. By the same reasoning, ()117p 
does not give a pinch surface, if 7^ 0. There are nevertheless two sets of 
solutions. First, there are infrared solutions where one line has zero momen- 
tum. 



Second, there are collinear solutions, where two of the lines are parallel to 
one of the outgoing external particles. 



The physical pictures associated with typical infrared and collinear pinch 
surfaces are shown in Figs. ITCT c) and (d), respectively. In each case, there is 
physically realizable propagation between vertices. 

Now we observe that even though solutions to the Landau equations like 
Eqs. ()118|) and ()119|) give pinch surfaces, they still do not necessarily produce 
mass dependence that is relevant to the leading power of Q, and hence are 
not necessarily leading regions. 



= 0, «! = ^2 = 0, 
(p + kf = 0, ai = as = 0, 
{p-kY = 0, 02 = 03 = 0. 




ai{p - ky + oi^ik^ = 0, 02 = 0, P = p ■ k = 0, 
Mp' + kY + asF = 0, ai = 0, k'^ =p' ■k = 0. 



(119) 
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The Bom graph for the vertex behaves hke Q". The contribution to the 
one-loop graph from the pure short- distance region, from the region, = 
0{Q), is Q"'~^. Thus this region is leading when the theory is renormalizable, 
at n = 6, but is non- leading relative to the Born graph when the theory is 
super- renormalizable, n < 6. 

Next we consider the one-loop graph near its singular surfaces. For exam- 
ple, consider the integral (jll5j) near the surface defined by the first of Eqs. 
()118p . To be specific, let j/c'^l < fcmax, /i = ■ ■ ■ 3, with fcmax being some fixed 
scale (which must > m). In this region the integral behaves as 



Compared to the short- distance region, this infrared region is leading only 
for n < 4. (The other two infrared regions in (jll8|) require n < 2.) (We 
remind the reader that, compared to the Born graph, none of these regions 
contribute to the leading power of Q.) 

Similarly, near the collinear pinch surfaces of Eq. ()118|) . the integral be- 
haves as 



so that again only for n < 4 do we find collinear contributions from this 
diagram that are leading compared with the short- distance contribution. 

In summary, for the scalar theory in six dimensions, only short distance 
regions are (relatively) leading for Fig. ^1 This result generalizes to all orders 
in the vertex correction for this theory 

The process of estimating the strength of a singularity is known as "power 
counting" . We will give more low order examples below, while more general 
arguments can be found in Ref. j2Sl- We can, however, summarize the ba- 
sic result of these arguments briefly. Let D be a reduced diagram with S 
"infrared loops" and Is infrared lines whose momenta vanish at the corre- 
sponding pinch surface, and with C "collinear loops" and Iq collinear lines 
whose momenta become proportional to an external momentum at the pinch 
surface. Finally, let denote the number of two point subdiagrams in R. 
In 0^ theory in n space-time dimensions, let us define the "degree of 
divergence" by 




(120) 




(121) 



uj{R) 



nS + {n/2)C - 2Is - Ic + N^. 



(122) 
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This corresponds to a power law 

Q{n-6)V ^<.{R) (^23) 

as Q ^ oo. Here V is the number of loops and A is a small parameter that 
parameterizes the approach to the singularities in the massless theory. The 
power law is measured relative to the power for the Born graph. 

In the renormalizable case, n = 6, we get leading behavior only if uj{R) = 
and there are no graphs that give negative uj{R). The term N2, which 
tends to suppress the behavior of the integral at the singular point, is due 
to the fact that the renormalized two-point function must vanish on-shell if 
the particle is to have zero mass. In this sense, the infrared behavior of the 
theory is dependent on renormalization 

In the superrenormalizable case, n < 6, the first factor in ()123|) gives a 
negative power of Q that just corresponds to normal ultraviolet power count- 
ing; this is the power that comes from a purely short-distance contribution to 
the graph. A sufficiently strong power law singularity in the massless theory 
is needed to overcome this if one is to get a leading contribution. 

7.5 Leading regions for deeply inelastic scattering in 

As an application, we now discuss the general leading regions for the basic 
inclusive cross sections in (0^)g. The criterion of physical propagation shows 
why the considerations of Sects. ElandEltake into account all relevant leading 
regions. To see this, we must generalize our concept of leading regions to 
include cut diagrams, of the type discussed in Sects. El and El There is no 
problem in doing this, and power counting may be estimated for collinear 
and infrared cut lines with the same degree of divergence Eq. (jl22j) as for 
virtual lines. 

It is useful to apply the optical theorem to reexpress the deeply inelastic 
scattering structure function Eq. as the discontinuity of the forward 

Compton scattering amplitude T{q,p), 

F{x,Q'^) = discT(g,p), 

nq,p) = % j d^l/e"'-^(p|T(j(y)j(0)b). (124) 
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Figure 17: (a) Diagrams for forward Compton scattering, (b) Reduced dia- 
grams for pinch surfaces, (c) Reduced diagram for a singular surface which 
is not pinched. 

This relation holds diagram-by-diagram, once cuts are summed over, so that 
a necessary condition for a region L to be leading in F is that it be leading 
in T. 

We thus need to consider the leading regions of the diagrams illustrated by 
Fig. lTTT a). which represents forward Compton scattering. The pinch surfaces 
are symbolized in Fig. ITTT b). The incoming hadron can form a jet of lines, 
and in addition, may interact with any number of soft lines, connected to a 
subdiagram 5* consisting of only lines with zero momentum. To see why this 
is the general form, consider a singular surface not of this kind, as in Fig. 
ITTT c). Here, one or more lines of the jet may scattering with the incoming 
photon to form a set of on-shell outgoing jets, which then rescatter to emit the 
outgoing photon and reform the outgoing jet, which eventually evolves into 
the outgoing hadron. Such a process is certainly consistent with momentum 



59 



conservation. Surfaces of this type are not pinch surfaces for the amphtude 
T{q,p), however, because the outgoing jets can never colhde again once they 
have gone a finite distance from the point at which they are produced. As 
a resuh, in every pinch surface, the incoming and outgoing photons attach 
at the same point in space-time, and we derive the picture of Fig. ITTT b). 
This resuh shows that all divergences associated with final state interactions 
cancel in the sum over final states. 

As indicated above, not every pinch surface will correspond to a leading 
region. In particular, power counting using Eq. ()122j) shows that there are 
no infrared divergences in {4>^)q, that is, no zero- momentum lines for any 
leading region [231 • In addition, we can show that only the minimum number 
of jet lines (two) can attach the jet to the hard part in Fig. ITTT b). It is a 
straightforward exercise in counting to prove these results, using ()122j) . the 
Euler identity (loops = lines - vertices + 1), and the observation that every 
internal line of a graph begins and ends at a vertex. In summary, we can 
show that Fig. 121 is indeed the reduced diagram of the most general leading 
region for deeply inelastic scattering in the scalar theory. 

7.6 Unitarity and jets: the cancellation of final state 
interactions 

The cancellation of final state interactions in deeply inelastic scattering plays 
an important role in the analysis for deeply inelastic scattering just described. 
This cancellation is a general feature of inclusive hard scattering cross sec- 
tions, and is used repeatedly in factorization proofs. The physics behind this 
cancellation has already been pointed out in Sect. El a hard scattering is well 
localized in space-time, and, as a result, it cannot interfere with long-distance 
effects which describe the further evolution of the system. Thus, when we 
sum over final states in an inclusive cross section, we lose information on the 
details of evolution in the final state, and are left with the constraint that, by 
unitarity, the sum of probabilities of all final states is unity. As a result, at 
each order in perturbation theory, long-distance contributions to final states 
must cancel. 

It is worth noting that it is not always necessary to sum over all final 
states to cancel long-distance interactions. There are three kinds of cross 
sections, among those mentioned in Sect. ^ for which the cancellation of 
final state interactions is important. In deeply inelastic scattering and Drell- 
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Yan, for instance, we sum over all hadronic final states. In single-particle 
inclusive cross sections, on the other hand, we shall find in Sect. IHl that 
cancellation requires the use of Ward identities. Finally, in jet cross sections, 
the cancellation comes about in a sum over all final states which satisfy 
certain criteria in phase space. Let us hint at how this happens. 

If all non-forward particles in the final state emerge from a single hard 
scattering, the criterion of physical propagation requires that the long-distance 
contributions will come entirely from soft and collinear interactions. This is 
because, as in the low order example of Fig. ^1 jets emerging from a sin- 
gle point and propagating freely cannot meet again to produce a new hard 
scattering. In this case, once the energy and direction of a set of jets is spec- 
ified, the sum over only those final states consistent with these jets will also 
give unity, and their collinear and infrared divergences cancel in the sum. 
The technical proof of these statements may be given in a number of ways. 
The simplest is based on a truncation of the hamiltonian to describe only 
collinear and infrared interactions. Then, since the truncated hamiltonian is 
hermitian, it generates a unitary evolution operator whose divergences cancel 
by the "KLN" theorem ^Hl 1^ . It is also useful to see that this cancellation 
is manifested on a diagram-by-diagram basis within each leading region in 
perturbation theory |23]- Proofs of this type are most easily given in terms 
of time-ordered or light-cone ordered perturbation theory [TUl I27j. 

Technicalities aside, the cancellation of final state interactions at the level 
of jets has a number of important consequences. The simplest of these is the 
finiteness of jet cross sections in e"'"e~-annihilation cross sections [2H1- We 
have already seen its importance for the analysis of deeply inelastic scattering 
in {(P^)q- It has a similar simplifying effect for single-particle inclusive cross 
sections, as well as for the Drell-Yan and related cross sections. To illustrate 
this, we show, in Fig. the reduced diagrams for leading regions in the 
scalar "Drell-Yan" cross section, defined for the scalar theory by analogy to 
Eq. dnSl) in (0^)6, after the sum over final states. We see that all information 
about the final state has been absorbed into a single hard part H. Note that 
this result holds not only for the fully inclusive Drell-Yan cross section, but 
also for semiinclusive cross sections such as hadron-hadron Drell-Yan pair 
-|- jets. 
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Figure 18: Leading regions for the Drell-Yan cross section in 

8 Factorization and Gauge Invar iance 



In this and the following section, we discuss the extension of factorization 
proofs to gauge theories. We begin with a discussion of the classical Coulomb 
field of a fast moving charge, an example that anticipates what happens in 
the full quantum theory. 

We next summarize the Ward identities we will need. Then we discuss 
the leading regions of Feynman graphs in a gauge theory. There are great 
differences from the case of the (p^ theory discussed in the previous section. 
With the aid of the example of the vertex graph, we show how, after an 
appropriate eikonal approximation. Ward identities are applicable that will 
combine graphs into a factorized form. 

In the next section, we will show how factorization may be proved for a 
variety of experimentally important cross sections, which can be measured 
in deeply inelastic scattering, e~^e~ annihilation and hadron-hadron scatter- 
ing. We should emphasize at the outset that although we regard existing 
proofs in all these cases as reasonably satisfying, there is still room for im- 
provement, especially for hadron-hadron cross sections. We will point out 
the shortcomings of existing arguments in Sect. ^1 



8.1 Classical considerations 

Before getting into a detailed discussion of Feynman diagrams, it is worth 
noting that insight can be gained into the physical content of factorization 
theorems from purely classical considerations. This discussion will at once 
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highlight an important difference between gauge and scalar theories, and at 
the same time show why this difference, important though it is, respects 
factorization. 

As we observed in Sect. 11.41 the parton model picture of hadron-hadron 
scattering rests in part on the Lorentz contraction of colliding hadrons. Now 
a simplified classical analog of a hadron is a collection of point charges, each 
acting as a source of a classical scalar field. We would expect that if the 
parton model, or factorization, is to make sense, these fields ought to be 
Lorentz contracted themselves, and this is just what happens. Let us see 
how. 

Consider first a static classical scalar field 0ci(x), associated with a point 
particle of charge q at the origin. If we assume that the field obeys Laplace's 
equation, it is given in the rest frame of the particle by 

0ci(x) = ^ (125) 

Now consider the same field in a frame where the particle is moving at velocity 
c/3 along the z-axis. Then the field at x'^ in this frame is 

0ci(a;') = T72^ (126) 

where, as usual, 7 = (1 — /3^)~^/^. For an observer at t' = in the primed 
system, the field decreases as 7"^ as P approaches unity, except near Xg = 0. 
Thus, the field is indeed Lorentz contracted, and any force proportional to 
the (j) field is also Lorentz contracted into a small longitudinal distance about 
x'^ = 0. This means that in the rest frame of a scalar "hadron", the forces 
due to another such hadron approaching at nearly the speed of light are 
experienced in a Lorentz-contracted fashion, just as supposed in the parton 
model. 

Now let us apply this reasoning to a classical gauge theory, in this case 
classical electrodynamics. Here, the field in the rest frame of a point particle 
of charge q is precisely analogous to Eq. 



<i(x) = (127) 

Because this is a vector field, however, there is a big difference from the scalar 
case in a frame in which the particle moves with velocity 0/^x3. In this frame. 
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we find 



KitW) = 0. (128) 

For large 7, the field in the zero and three directions are actually independent 
of 7 at fixed times before the coUision. It might therefore seem that a vector 
field is not Lorentz contracted, and would not respect the assumptions of 
the parton model. If we look, however, at the field strengths rather than 
the vector potential, we find a different story. The electric field in the three 
direction, for instance, is given in the primed frame by 



^aM = . r'r: (129) 



which shows a 7"^ falloff. Since the force experienced by a test charge (or 
parton) in the primed frame is proportional to the field strength rather than 
the vector potential itself, the physical effects of the moving charge are much 
smaller than its vector potential at any fixed time before the collision. This 
in turn may be understood as the fact that, as 7 00, the vector potential 
approaches the total derivative 

qd^'lnipct' -x's). (130) 

That is, for any fixed time the vector potential becomes gauge equivalent to 
a zero potential. 

We can conclude from this excursion into special relativity that factor- 
ization will be a more complicated issue in gauge theories than in scalar 
theories. Only for gauge invariant quantities will the gauge-dependent, large 
vector potentials of moving charges, which naively break factorization, cancel. 
So, in particular, we cannot expect factorization to be a property of individ- 
ual Feynman diagrams, as it was in scalar theories. On the other hand, we 
should look for the solution to these problems in the same techniques which 
are used to show the gauge independence of physical quantities. 
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Figure 19: Ward identity. 
8.2 Ward identities 

The Ward-Takahashi identities of QED and the Taylor- Slavnov identities of 
nonabehan gauge theories ensure the perturbative unitarity of these theories. 
We shall refer to them collectively as 'Ward identities" below. 

Ward identities may be expressed in various forms, for instance, as iden- 
tities between renormalization constants (the familiar Z\ = Z2 of QED). For 
our purposes, however, the basic Ward identity is given graphically by the 
equation 

{N\ T d^,A^'{xi) X ■ ■ ■ X |M) = 0, (131) 

where A'^{x) is an abelian or nonabelian gauge field, and where M and are 
physical states, that is, states involving on-shell fermions and gauge particles, 
all with physical polarizations. In particular, physical states do not include 
ghosts. Equation ()131|) will be represented graphically by Fig. ^[ in which 
the scalar operator d^A'^lx) is represented by a dashed line ending in an 
arrow. In momentum space, this operator is associated with a standard 
perturbation theory vertex in which one gluon field is contracted into its 
own momentum. Here and below, we refer to such a gluon as "longitudinally 
polarized" . Note that this is to be taken as referring to the four-momentum. 

Proofs of Eq. ()13H) are most easily given in a path integral formulation 
using BRST invariance, as in, for instance Ref . [2^] . They can also be proved 
in a purely graphical form, as in the original proofs of Refs. jHHI and p?T] . 
Here we need not concern ourselves with the details of these proofs, although 
it may be worthwhile to exhibit the very simplest example of Eq. ()131|) . This 
is the lowest order contribution to the electron scattering amplitude with a 
single longitudinally polarized photon. At this order, we have 

q'^u{p + q)'^^u{p) = u{p + q)[{^ + + m) - {■/) + ■m)]u{p) = 0. (132) 

The first equality is sometimes referred to as the "Feynman identity" , and the 
overall result is current conservation at lowest order. This is not surprising, 
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since classical current conservation is a consequence of gauge invariance. In 
the quantum theory, it appears as a matrix element relation, whose validity 
is ensured by the Ward identity. A helpful exercise is to construct the analog 
of Eq. ()132|) for the scattering of a physically polarized gluon. The graphical 
proof consists essentially of repeated applications of identities like Eq. ()132|) . 

Even without going into the details of the proof of Eq. (jl31|) . we can 
elucidate its interpretation. First, it is true order-by-order in perturbation 
theory, although not graph-by-graph in perturbation theory. In addition, we 
may imagine constructing a path integral in which only certain momenta 
are included, for instance ultraviolet momenta and/or momenta parallel to 
a given direction. Then at a given order, the Ward identities hold for both 
internal and external lines in this restricted portion of momentum space. 
This heuristic argument may be verified by a close look at the graphical 
proof of Ward identities in Refs. jHO] and [SI] . 

So far, we have discussed Ward identities for the S-matrix. As we saw 
in Sect. lU however, we will sometimes be interested in matrix elements in- 
volving a gauge invariant but nonlocal operator which includes the ordered 
exponential of the gauge field. Such matrix elements also obey Ward identi- 
ties, which may be proved by either of the methods mentioned in connection 
with Eq. ()131|) . The simplest generalization of ()131|) to this case is 



{N\T {U,d^A>^'{x,) Vexp{zg / d2/-A+(0, y", Ox)}$(0)) |M) = 0, (133) 



where in the ordered exponential A refers to the gauge field in the represen- 
tation of field which may be a fermion or gauge particle. Equation ()133|) . 
in various guises, will be useful in our proofs of factorization. 

8.3 Singularities in Gauge Theories 

Discussions of factorization start with a catalog of the pinch surfaces of the 
relevant Feynman diagrams, as described in Sect. [7| They then proceed, by 
power counting, to estimate the strength of singularities encountered in each 
such surface. The same procedure may be carried out for gauge theories, 
but, as we will now see, many of the regions that are nonleading in (p^ are 
now leading. Thus the results are much richer than in 0^ theory. 

The one- loop vertex graph illustrates the origin of the infrared and collinear 
singularities. Ignoring overall factors, including group structure, we find that 
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the graph is given by 




(27r)4 Kp' - ky - m2 + ie] [{p + kf - w? ^ %e\ {k^ + ie) ' 



(134) 



The singularity structure (with one minor exception) is the same as in (j)^ 
theory; what changes is the strength of the singularities. To discuss the large 
Q region, we will consider, as before, the massless limit. Of the solutions 
(I118|) and (jll9|) to the Landau equations, the first of the infrared solutions 
{k^ = 0) and both of the coUinear solutions (jll9p give leading power behav- 
ior at large Q, as we will see. These are singularities in the fully massless 
theory, and, by our discussion in Sect. [7J they correspond to long-distance 
contributions when Q is large. 

In addition to these singularities, there is a genuine singularity, at k'^ = 0, 
even when the fermion mass is nonzero. This is an example of the usual in- 
frared divergences of QED and is caused by the masslessness of the gluon. 
This singularity survives the Q —>■ oo limit, of course, and becomes the first 
of the solutions (jllSj) . The methods that we use to treat both the coUinear 
and especially the infrared singularities in the fully massless nonabelian the- 
ory are explicitly motivated by the elegant methods given by Grammer and 
Yennie jS21 to treat the ordinary infrared problem in QED. In QED, the in- 
frared divergences correspond directly to the real physics of the long range 
of the Coulomb field and the genuinely massless photon. But in QCD the 
infrared divergences are cut off by confinement. Since this is a nonpertur- 
bative phenomenon, the resulting cutoff is not easily accessible (if at all) in 
perturbation theory. Perturbative calculations must be restricted to suffi- 
ciently short-distance phenomena so that asymptotic freedom is useful. The 
singularity structure of the massless theory is just a convenient tool to aid 
in the factorization of long-distance phenomena. 

8.4 Infrared divergences in gauge theories: the eikonal 
approximation 

We now consider the infrared singularity at A;^ = 0. Although our ultimate 
aim is to treat the large Q limit, our discussion will not need to assume this 
limit initially. Let us see how the integral ()134j) behaves near this point. As 
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— s> 0, it is valid to make the following two approximations in V^(p, p'); 

(1) Neglect k^^ compared to m and p'^ in numerator factors, 

(2) Neglect fc^ compared to p' ■ k and p' ■ k in denominator factors. (135) 

Together, these two prescriptions define the "eikonal approximation" for the 
graph. Simple manipulations show that in the eikonal approximation is 
given by 

V,ip,p') = -Ap.p'vip'hAp) / i-2Jf.k + ^e){2lk + ^e){k^ + ^ey 

(136) 

In this form it is apparent that the k integral is logarithmically divergent 
from the region near k^ = 0. Notice also that, because the numerator in 
proportional to p' ■ p, this divergent integral behaves as a constant at as 
p' -p oo, that is, with the same power as the elementary vertex. This is to 
be contrasted with the situation in 0^ theory that was explained in Sect. [7| 
Infrared behavior with the same power law behavior in Q as the elementary 
vertex is a characteristic of theories with vector particles. 

The eikonal approximation is, not surprisingly, closely related to the or- 
dered exponentials of Sect. 14. 2t with their eikonal Feynman rules. In fact in 
making the eikonal approximation p35|) . we are precisely replacing fermion 
propagators by eikonal propagators of the type shown in Eq. (P5|) for the par- 
ton distribution functions. We can anticipate the importance of the eikonal 
approximation by relating it to the classical discussion given in Sect. 18.11 
Consider a gluon of momentum fc^ interacting with an eikonal line in the 
y/j. direction. The only component of the gluon momentum k^ on which the 
eikonal propagator depends is v ■ k, and the only component of the gluon 
polarization e'^{k) on which the eikonal vertices depends is f ■ e. So, as far 
as the eikonal line is concerned, the gluon acts in the same way as a ficti- 
tious gluon of momentum {v ■ k)u^ and polarization {v ■ e)u^, where is 
any vector for which u ■ v = 1. But this fictitious gluon is longitudinally 
polarized. That is, any gluon interacts with an eikonal line in the same way 
as a longitudinally polarized, and therefore unphysical, gluon. But we have 
argued above that such gluons, although they may be expected to break fac- 
torization on a graph-by-graph basis, should be consistent with it in gauge 
invariant quantities. 

When Q is large, we consider not just the actual infrared singularity at 
k^ = 0, but the whole infrared region fc^ <^ Q. That is, as Q — > oo, we 
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consider the region k^/Q 0. It is possible for the different components 
of k^^/Q to go to zero at such different rates that the eikonal approximation 
fails. Since we will rely on this approximation in proving factorization we 
will need to evade this failure. 

To get an idea of what is involved, let us return to Eq. ()134p . and justify 
the eikonal approximation Eq. (jl35p in this simplest of cases. Failure of Eq. 
(jl35|) is caused by failure of the second of the approximations of which it is 
comprised: dropping factors of k^ in the numerator is a safe bet, because, as 
we have seen, the factors combine to form large invariants. So, the issue 
is whether or not we may neglect k"^ compared to p ■ k and p' ■ k. This is 
nontrivial, because it is easy to find vectors k^ for which p ■ k and p' ■ k are 
small, while k'^ remains relatively large. This will be the case whenever its 
spatial momentum transverse to the p and p' directions is large, 

k^^ k'^ r^p-k, p' -k. (137) 

This region was called the "Glauber" region in Ref. jnS|- It is easy to check 
that in this region the k^ integral of Eq. (|134p is logarithmically divergent. 
If we were to put in a gluon mass (as is consistent for an abelian theory), the 
divergence would disappear, but we would still have a contribution from the 
region ()137p to the leading-power behavior at large Q, that is, a contribution 
of order times logarithms. 

Does this mean that the eikonal approximation is wrong? In fact it does 
not. To see this, we appeal to our freedom to deform momentum space 
contours. Suppose, for definiteness that p and p' are in the ±z directions, 
respectively, and that |p| = |p'|. We then change variables from the set 
{/co,kj} to the set {k^ = 2~^^'^(kQ ± (p ■ \i)/ujp), k±}. (These become light 
cone variables in the high energy limit.) Then in the region defined by Eq. 
(I137p . the fc^ ^ singularity of (|134p is given by 

1 1 

d^/c I T -dn^dK' 



-kj_^ + ie) - k^^ + ie) (-23/2cu;k+ - k^^ + ie) ' 

(138) 

where the variables appear in only one denominator each. In this form we 
see explicitly that the integrals are not trapped in the region <C /cj_, 
since they each encounter only a single pole in this region. As a result, these 
contours may be deformed away from the origin into the region ~ 
But in this region the eikonal approximation is valid, provided only that 
\k^\ -C Q. So, we may relax our criteria for the eikonal approximation 
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to include the possibility that, even if it is not valid everywhere along the 
undeformed contours, these integrals can be deformed in such a way that it 
holds along the deformed contours. 



8.5 CoUinear divergences and choice of gauge 

In addition to infrared divergences, we have to consider coUinear divergences 
in the massless limit. The nature of the coUinear contributions to leading 
regions depends on the gauge, as we will now show. 

Consider the gluon propagator in a axial gauge n - A = 0. It has the form 

D^Ak) ^ ^ - ^ ^) , (139) 

which satisfies 

. . (i!^ _ ^) . (140) 

Such a gauge is "physical" because its propagator has no particle pole when 
contracted into any vector proportional to its momentum. Another way of 
putting this is that in such a gauge longitudinal degrees of freedom do not 
propagate. This is to be contrasted with a covariant gauge like the Feynman 
gauge, for which k^D^y{k) = —ik^/k'^. As a result, leading regions in which 
longitudinal degrees of freedom propagate are present in covariant gauges 
but absent in physical gauges. Let us see what this means in practice. To 
do so, we turn again to the vertex correction, (jl34|) . 

We have already stated the locations of the collinear singularities, in Eq. 
flll9p . The two possibilities are that k^ is proportional to and that it is 
proportional to p'^ . The corresponding reduced diagrams are shown in Fig. 
llOfa). 

By doing the k~ integral by contour integration, we easily find that in 
Feynman gauge the contribution of momenta close to the singularity where 
k is proportional to p is given by 

VMp') ^ y -^v{p)iMp) 9-^ (1 + (141) 

and similarly in the second region. We have exhibited explicitly the numera- 
tor of the gluon propagator. For fixed /c^ = —xp~^, the k± integral diverges, 
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(b) 



Figure 20: Leading coUinear reduced diagrams at one loop: (a) covariant 
gauge, (b) physical gauge. 

and is leading power, that is, independent of Q. This is the coUinear diver- 
gence. (There is an additional infrared divergence as fc^ vanishes; this region 
we have already discussed. This result for the Q dependence is known as 
a "Sudakov" double logarithm; it is associated with the overlap of coUinear 
and infrared divergences.) 

These regions, summarized by the reduced graphs of Fig. EUT a) in which 
two coUinear lines attach to a hard subdiagram, would not be leading in 
(0'^)g, because of the lack of the numerator factor. Note, however, that in the 
numerator of the gluon propagator, the term which gives the leading behavior 
in the collinear region is (?_+. Since the gluon is moving, by assumption, 
parallel to p^, which is in the plus direction, this corresponds to an unphysical 
polarization at the vertex adjacent to the antiquark line. Thus the collinear 
divergence is associated with a longitudinally polarized gluon, and we might 
expect it to be absent in a physical gauge — at least in this particular 
diagram. 

To verify this, we can compute in an axial gauge. The leading term 
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in Eq. (|14ip is then replaced by 



/d'^k± / n^k+ + A;_A;-|_ 

Using Eq. ()14()j) (and remembering that k± = k^), we easily check that the 
collinear divergence in ()141|) is absent in ()142|) . and that the vertex diagram 
therefore lacks the Sudakov double logarithm in axial gauge. Of course, since 
the theory is gauge invariant, the corresponding physics, and in particular 
the double logarithm, has to show up somewhere, and in axial gauge it occurs 
in the one loop fermion self energy. We leave it as a simple but instructive 
exercise to check that this is indeed the case. Thus, in axial gauge the reduced 
diagrams of Fig. I^UT a) do not correspond to a collinear divergence, while that 
of Fig. EUT b) does. We emphasize here the fact that in the Feynman gauge the 
jets are one-particle irreducible, while in the axial gauge they are reducible. 
As we shall see below, this result generalizes to all orders. In this sense, the 
gauge theory in a physical gauge behaves, from the point of view of reduced 
diagrams for collinear lines, like (0^)g. 

This suggests that an axial gauge is the most appropriate one for proving 
factorization. However, the singularities at n ■ A; = cause a lot of trouble. 
In the first place they obstruct ESI ISn] the contour deformations that 
we have already seen are essential to demonstrating factorization; this is 
equivalent to saying that the singularities violate relativistic causality on a 
graph-by-graph basis. Furthermore, the singularities have to be defined by 
some kind of principal value prescription, and it is difficult to ensure that 
the products of these singularities that occur in higher order graphs can be 
defined properly |?7] . 

8.6 Power counting for gauge theories 

As in the scalar theory, we must use power counting to identify those pinch 
surfaces which actually give leading regions. Again, this approach is discussed 
in detail in Ref. [^. Here, we once again quote the general result. Assuming 
the eikonal approximation, for any leading region with reduced diagram R, 
we compute the infrared degree of divergence, uj{R), analogous to Eq. ()122|) . 
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coiR) = AS + 2C- 2Is -IC + N2 + ^Ns + F. (143) 

We have assumed a space-time dimension 4. As in ()122p . 5, Is, C, and Ic 
are, respectively the numbers of soft loops and lines, and coUinear loops and 
lines in R at the associated pinch surface. N2 is the number of two-point 
functions in R, while N3 the number of three-point functions all of whose 
external lines are in the same jet. F is derived from the numerator factors 
where soft lines connect to coUinear lines. It is positive except when all soft 
lines connecting to coUinear lines are gluons. The suppression terms, ^N^ 
and F, are the only differences from Eq. (|122|) . 

We note that the term is present diagram-by-diagram in physical 
gauges 13 , but that in Feynman gauge the computation holds in general only 
when gauge invariant sets of diagrams are combined for the hard scattering 
subdiagrams. In fact, in covariant gauges, individual diagrams may be much 
more divergent in the presence of infrared and coUinear interactions than is 
the cross section, and may even grow with energy jHSj- This is a consequence 
of the well-known fact that unitarity bounds on energy growth are only a 
property of gauge invariant sets of diagrams. 

8.7 General leading regions 

The general leading regions for e"''e~-annihilation processes, for deeply in- 
elastic scattering and for hard inclusive hadron-hadron scattering are quite 
analogous to those for the (0^)g theory. The basic difference is that lines 
which participate in infrared logarithms must be added to the corresponding 
reduced diagrams. 

Fig- I^Tl a) shows a general leading region for a single particle inclusive 
cross section in e~^e~ annihilation for a physical gauge, and Fig. l?1T b) for a 
covariant gauge. 

Compared to the leading regions for (0^)g, summarized in Fig. which 
include only jets of coUinear lines and hard subdiagrams, to get Fig. l^lT a) we 
simply add a "soft" subdiagram, consisting of lines whose momenta vanish 
at the pinch singular surface in question. The soft subdiagram contains in 
general both soft gluon lines and soft quark loops (as well as ghost loops 
in covariant gauges); its external lines, however, are always gluons. These 
external gluons always attach to (energetic) lines and not to the hard subdia- 
grams. The lines attaching jet subdiagrams to the hard subdiagrams may be 
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Figure 21: Typical leading regions for annihilation processes, (a) physi- 
cal gauge, (b) covariant gauge. The most general leading region has the 
possibility of extra jets beyond the two shown here. 
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either gluons or fermions, but at leading power only a single line from each jet 
enters a given hard subdiagram, just as in (0^)g. The physical picture is also 
the same as in {4>^)q', several hard particles recede from a hard scattering at 
the speed of light, forming jets by their self-interactions. These particles can 
never interact with each other except by transfer of soft momenta \k^\ <^Q. 
The presence of vector particles in the gauge theory, however, does give lead- 
ing power contributions from the exchange of soft particles. Since /Q ~ 
for each of these particles, they do not affect those of the Landau equations, 
flll2p and ()113|) . which involve only the jet subdiagrams. Of course, it should 
be kept in mind that the soft lines have zero momentum only at the exact 
pinch singular surface. Feynman integrals get contributions from an entire 
region near this surface where the soft momenta are much smaller than a 
typical energy of a jet, but may approach a nonzero fraction of that energy. 

For a covariant gauge, the leading regions are essentially the same, ex- 
cept that, just as in the one-loop case, arbitrary numbers of longitudinally 
polarized gluons may attach the jet to the hard part, as shown in Fig. [?lT b). 

Figures YFlXa) and E^b) show general leading regions for inclusive deeply 
inelastic scattering, and the Drell-Yan cross sections in Feynman gauge. As 
with (0^)g, the sum over final states eliminates pinch singularities involving 
final state jets. The remaining on-shell lines make up the jets associated with 
the incoming particles including soft exchanges within and between the jets. 
In covariant gauges, longitudinally polarized gluons may connect the jets to 
the hard part. 

9 Factorization Proofs in Gauge Theories 

We are now ready to discuss the extension of factorization theorems to gauge 
theories, for the basic cross sections discussed above: deeply inelastic scat- 
tering, single particle inclusive annihilation and Drell-Yan production. Each 
of these will require new reasoning relative to the scalar case. 

Compared to the proof in Sect. Elfor (0^)g, our treatment of factorization 
in gauge theories will be much more modest. Rather than derive closed 
expressions for the factorized forms in terms of explicit subtraction operators, 
we will deal with the cross sections on a region-by-region basis. We will show 
that an arbitrary leading region either contributes to the factorized form of 
the cross section, or cancels to leading power when gauge invariant sets of 
diagrams are combined. 



75 




Figure 22: Leading regions in Feynman gauge for (a) inclusive deeply inelastic 
scattering and (b) Drell-Yan cross sections. 
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9.1 Deeply inelastic scattering and collinear factoriza- 
tion 

We start with the deeply inelastic scattering cross section, h{p) + (q) ^, 
with p^{q'^) being the momentum of the incoming hadron h (virtual photon 
7*). Here, as we shall see, the question of factorization reduces to a treatment 
of collinear singularities associated with unphysically polarized gluons. 

Fig- I22ta) illustrates the leading regions in a gauge theory for diagrams 
that contribute to the structure function tensor W^'^{q,p). There is a single 
jet J, in the direction of the incoming particle, and a single hard subdiagram 
H'^'^, containing the hard scattering. Divergences associated with final state 
interactions cancel because of unitarity in the sum over different final state 
cuts of the same Feynman graph. Thus we have not included regions in which 
soft gluons from the jet J interact with the outgoing particles in the hard 
part H, even though such regions can give leading contributions to individual 
cut graphs. 

In Feynman gauge, the hard subdiagram is connected to the single jet by 
more than one collinear line. This makes the transition to the convolution 
form of Eq. (j2)) more complex than in the scalar case. 

In physical gauges, the reduced diagram corresponding to an arbitrary 
leading region has the same form as for the scalar theory. Fig. El This sim- 
plification is the reason that most of the original arguments for factorization 
were given in physical gauges 14, where essentially the same procedure can 
be used as in Sect. However, it is important to show how the proof may 
be carried out in the covariant gauges, for two reasons. First, as mentioned 
above, there are difficulties associated with the unphysical singularities en- 
countered in physical gauges, which have not been fully understood yet . 
Although these are presumably of a technical nature, and not associated 
with the content of factorization, it is surely desirable not to be completely 
dependent on this presumption. Second, physical gauges, because of their 
noncovariance, are ill-suited to proofs of factorization in the crucial case of 
hadron-hadron scattering. So, in the interests of generality, we shall discuss 
deeply inelastic scattering in the Feynman gauge. These issues were not 
treated in Ref. 0. 

Let us consider a typical cut Feynman diagram G^'"^ where C labels the 
cut, in the neighborhood of a leading region L. L is specified completely once 
we specify how the graph G is to be decomposed into the subgraphs J and 
H. We shall write the contribution from region L to G^'"^ as G^^'^\ 
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Referring to Fig. E^ a). we see that our problem is to organize the set 
of longitudinally polarized lines which attach the jets to the hard parts. 
Suppose a set of n gluons of momentum /f' attaches to the hard part H to 
the left of the cut, along with a physically polarized parton of momentum 

— /f . Similarly, suppose a set of n' longitudinally polarized gluons 
/'^^ attaches to H on the right of the cut, along with a physically polarized 
parton of momentum k^ — 'Y^. I'j. Each momentum is parallel to the 
external momentum p'^ and flows into the hard part. Each I'j is also parallel 
to p^, but flows out. We sum over all cuts of the original graph G consistent 
with this leading region, with fixed n and n'. We can now represent the sum 
over these allowed cuts C of G^^'^^ as 

c 

= I ^ n /^n ! ^ 

J (27r)4 1.1 7 (27r)4lly (27r)4 

^ 3 

X J2 J^''-'\P''-^ - {r }; - Sr;. {^t (144) 

Cj 

where //j and Vj are polarization indices for the U and /'j, respectively, and 
Tj and f]' are the polarization indices associated with the physical partons 
attaching to the hard part on either side of the cut. Of necessity, the sum 
includes only those cuts which preserve n and ra', and we note that it breaks 
up into independent sums over the cuts of the hard part and of the jet. 

The integrals in (jl44j) are restricted to the neighborhood of the region L. 
We implicitly introduce a variable /x, to set the scale of L. The integration 
region in ()144j) is set by requiring, for instance, that lines within H have 
transverse momenta of order at least /i, while those in J have transverse 
momenta of /i or less. yU will later be identified with the renormalization 
scale for the parton distribution. 

Because all lines in H are, by construction, far off the mass shell, we 
replace the momenta of all its external particles by lightlike momenta in 
the corresponding jet direction. Then, if we keep only leading polarization 
components, the extra collinear gluons which attach the hard part to the jet 
are exactly longitudinally polarized. Corrections are suppressed by a power 
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of q^. To formalize this approximation, we introduce the vectors 



= gl, = gt, (145) 

and define 

u-11 = \i, u-k'' = k, u- C = (146) 
In terms of these variables, the approximation is 

Ch 

^ ^(g^ {k - SA,)t;^ {A,^;"»}; (A; - SA;.)^;^ {A^^ }),,,' J]"''' H^'^^ 

« i 

(147) 

where 

H{q'^- {k - SA,)t;^ {A,t;°»}; (k - SA;.)t;^ {A^^^^),,,. 
= J2 H^^"\q''\ kv^ - SA.^;^ {A,t;"^}; A;t;^ - SA;^;^ {A>^^})J;;;'^^ 

X 11% 11% 

This replacement is analogous to the operator P introduced for the scalar 
theory in Sect. IHl 

We will now show that the unphysical polarizations of the extra gluons 
can be used to factor them from the hard part. The hard part will become a 
function of only the total longitudinal momentum flowing between it and the 
jet, as is appropriate for a factorized form, while the longitudinally polarized 
gluons will couple to an eikonal line, which we associate with the jets. 

Let us show this result first for the diagram on the left-hand side of 
Fig- ESfa), with a single longitudinally polarized gluon of momentum Z'^, 
which attaches to the hard part along with the physically polarized parton of 
momentum k^ — l^. (We shall refer to particles by their momentum labels.) 
If we apply the Ward identity, Eq. p3H) to this set of diagrams, we find 
the result on the right hand side of Fig. |22ta). The left-hand side of Fig. 
I2nta) would vanish, except that the diagram on the right-hand side, in which 
the gluon is attached to the physical parton, is not included in H by 
construction. But now consider the identity shown in Fig. |22Ib). Here we 
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Figure 23: Ward identities for a single gluon. Group sums follow repeated 
indices, (a) Identity for hard part; (b) Eikonal identity; (c) Factorization of 
the gluon. 

consider a diagram in which the unphysical line ends in an eikonal line, while 
H has a single (physically polarized) external line from J, which carries the 
total jet momentum k. The right hand side of Figs. EHT a) and EHT b) are the 
same, and we derive the identity of Fig. EHT c). in which the longitudinally 
polarized gluon has been factored onto an eikonal line moving in the opposite 
direction from the A-jet. 

To be careful, we should note that in each individual cut diagram of Fig. 
1221 the intermediate states are not physical states, but rather states including 
on shell gluons with unphysical polarizations and ghosts. Once graphs for a 
given cut are summed over, however, we may replace the unphysical states 
by physical ones [30j. So we may, without loss of generality, treat the matrix 
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elements as though they were between physical states. 

The extension of this reasoning to two gluons is straightforward. We use 
the identity of Fig. |23 analogous to Fig. 1221 On the right hand side of the 
first equality in Fig. we have two diagrams in which only the physical 
parton attaches to the hard part, and also two diagrams in which one gluon 
is still attached to the hard part. (Diagrams in which the two gluon lines 
are interchanged are not indicated explicitly in the figure.) In a covariant 
gauge, Lorentz invariance requires that the gluon entering the hard part in 
diagram 4 also be longitudinally polarized (it has no other vectors on which 
to depend). Thus we can apply the result of Fig. 1221 for the single gluon 
entering the hard part in diagrams 3 and 4. The result is shown in the 
second and third equalities in Fig. (211 This inductive approach can clearly 
be extended to arbitrary order, and we derive Fig. |2S1 for a general leading 
region. 

This gives the overall replacement 

Hiq-"- {k - SA,)i;^ {Ait;"'}; {k - SA;.)^;^ {A>^^}),,,' H ^''^ H"'^' 

i j 

X j(^-^)(p- p - s/f, {/-}; k^ - sr;, {/t})?,t.} 

Hiq", kv^\^, £{u, {A,})^^'>^*(n, {A;.})^^^> 

X j(^^)(p- p - s/f , {/-}; k^ - sr;, (149) 

where S{u^ is a lightlike eikonal line in the direction, coupled to 

n gluons Zj, and similarly S*{u^ {■^j})^'^^^ is an eikonal in the same direction 
coupled to the n' gluons I'y It is natural to group the eikonal lines with the 
jet, and to define (compare Eq. (f73|l ) 

^(0- ^ /|^j*«p..A.„-i,n/||n/|;i 

Cj 

The function J is linked to the remaining hard part H{q'^ , kv^) through 
only the variable ^ and the physical polarization indices rj and rj' . Using 
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Figure 24: Application of Ward identities to two coUinear gluons. 
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Figure 25: Factorization of coUinear gluons. 
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(fTHIH) and (ITin|) in (fMl . we have 



^^^n^) = I^^H{q'',^p-uvnrj,r,'J{0''''\ (151) 

where a; = —q^/2p ■ g as in Sect. H Thus in each leading region the cross 
section factorizes into an ultraviolet contribution times a contribution to the 
distribution of the physical parton which remains attached to the hard part. 
It is clear that we get every leading region for the parton distribution in this 
way. Note that when we sum over all leading regions, the perturbative sums 
for the hard part H and the factorized jet J are completely independent. 

Equation (jl5ip contains most of the physics of factorization for deeply 
inelastic scattering, but a few more steps are required to obtain the result 
(121). First, one argues using Lorentz invariance that for unpolarized incoming 
hadrons the hard part are both diagonal in the spin indices i], rj' . Thus we 
can sum over the spin of the partons leaving the jet part and average over 
the spin of the parton entering the hard part. This decouples the two factors 
in spin space. We then sum over all graphs and over the leading regions L 
for each graph. The result is 

W^^'~E ['^^a/A{^^fi)K''iq'',^P-uv^,fi). (152) 

^ Jx ^ 

Here a sum over parton types a is indicated, is the jet part summed over 
graphs, leading regions, and spins, and Ti. is the hard part summed over 
graphs and leading regions and averaged over the incoming spins. Both 
and TC depend on the parameter fi that sets the scale for the leading regions. 

We can relate the functions JF to the MS parton distributions fa/A in 
the following manner. We note first that we can carry out exactly the same 
factorization procedure for parton distributions defined as in eqs. (j43|) and 
(I44p as for the deeply inelastic scattering structure functions above. Then, 
in place of ()152p . we find 

fb/Aiv,f^) = ^ [ y-^ba(WO •^aM(^,/i), (153) 

where gba is some new hard part (a matrix in the space of parton types), 
while J^a/A is the same jet function as in ()152|) . Here we define the scale of 
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the leading regions to be the same as the renormahzation scale in the parton 
distribution, and we use the same notation for both. 

Using eqs. ()152j) and ()153|) . we find the desired result for the structure 
functions, 

W^''{q^,pn ~ E - M^^) K\q',vp^,fi,a,{fi)) . (154) 
where the hard part H^'^ is defined by the relation 

j2 [^^ Hr{q^,r]p'^,fi,asifi)) ^.a(wo = ^^(?^ " z^)- (155) 
b ^ 

It should also be possible to demonstrate this factorization in the more 
careful manner outlined for the scalar theory in Sect. IHl From ()15H) . the 
leading region L may now be represented by Fig. El the canonical form for 
deeply inelastic scattering found in the scalar theory and in physical gauges. 
Since the same construction may be carried out for any leading region, one 
could define a subtraction procedure for gauge theories analogous to the one 
for scalar theories. The subtraction operator for a leading region L then 
makes the replacement p47p for the hard part for the region. 

9.2 Single-particle inclusive cross sections and the soft 
approximation 

The leading regions for a single-particle inclusive cross section, e"*" + e~ — *• 
A{p) + X, were shown in Fig. EH There is a jet subdiagram J that describes 
the jet in which hadron A is observed. The hard subdiagram H contains two 
short distance interactions (one on each side of the final state cut) involving 
highly virtual particles, from which one or more jets of interacting collinear 
particles emerge. Once again, there are extra longitudinally polarized glu- 
ons connecting the jet J to the hard subdiagram H. More importantly, in 
contrast with deeply inelastic scattering, there is a soft subdiagram S that 
connects J to H. As a result, the factorization property fails on a graph-by- 
graph basis. 

We recall that in any given cut Feynman graph for deeply inelastic scat- 
tering there could be soft partons connecting to the hard subdiagram, rep- 
resenting soft interactions between on-shell particles as they enter the final 



85 



state. However, we argued that any leading region containing such soft in- 
teractions gives a cancelhng contribution when one sums over the possible 
final state cuts for a given Feynman graph. Unfortunately, this rather trivial 
cancellation mechanism does not work for single particle inclusive produc- 
tion The reason is that we are observing a particle in the final state 
rather than summing freely over all final states. Nevertheless, our aim will 
be to show that any leading region with a soft subdiagram connecting J to H 
cancels. The only remaining leading regions are analogous to those already 
encountered in the case of deeply inelastic scattering, in Sect. 19. H so that 
the proof of factorization sketched there carries over. 

We consider the contribution from a leading region L to a cut Feynman 
diagram. Each such cut diagram is decomposed into subdiagrams J, H, and 
5*. We now sum over all cut graphs containing the same number of lines 
connecting the parts J, H, and S on each side of the cut and call the result 
G. Our object is to show that, after summing in addition over where the 
lines go relative to the final state cut, ^ G can be rewritten in a factorized 
form in the high energy limit. 

In order to write the kinematic approximations, we pick lightlike vectors 
= in the p'^ direction and u'^ = gt in the Q^,p^ plane, and define the 
momentum fraction ^ of the outgoing hadron A by 

u ■ k = u ■ p/ ^. (156) 

As in the case of deeply inelastic scattering, we use the longitudinal po- 
larization of the extra collinear gluon lines which attach the J to the hard 
part H. We once again approximate these lines by dominant momentum 
and polarization components, so that they appear as longitudinally polar- 
ized. Then, we sum over graphs representing different attachments of the 
collinear gluons to H and use Ward identities to remove them from H and 
attach them instead to eikonal lines £ in the direction. They are then 
grouped with the jet to form, in this fragmentation function dA/a{C}- 

We thus derive a form analogous to ()151|) . but with the extra complication 
of the soft lines, 

xiJ(g^(n•p/Ot;^{f})''''^^^^ (157) 
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Here J is analogous to the function J in Eq. (|15(J|) . It includes the original 
jet subgraph together with eikonal lines attached to the 'extra' longitudinally 
polarized gluons that formally attached to the hard part. The indices 77,77' 
represent the polarization of the physical parton that enters the hard part 
carrying momentum A;^ ~ (u ■ p/S,)v^. The extra complication compared to 
deeply inelastic scattering is the soft subgraph S, which couples to J and H 
via soft gluons as indicated in (jl57|) . 

It is useful to interpret ()157|) in the language of Sect. 19.11 It describes the 
jet A containing the observed hadron A together with the unobserved jets 
that we have included in H. All of these jets emerge from a hard scattering 
and evolve independently, except for the exchange of soft partons, which are 
coupled to the color current of each jet. In the frame of jet A, for instance, 
all the charges within the other jets are moving at nearly the speed of light. 
But then, according to the discussion of Sect. 19.11 the Lorentz transformed 
field due to these jets should be nearly gauge equivalent to zero. Of course, 
since jet A arises from a quark or gluon, it is not gauge invariant. We might 
expect, however, that we can exhibit the gauge nature of this interaction. 
To do so, we need a generalization of the eikonal approximation which we 
applied in Sect. 18.41 to single parton lines coupled to soft radiation. 

The relevant generalization of the eikonal approximation has been termed 
the "soft approximation" |35| IH^ 1^ l4Uj . For the A-jet,it consists of making 
the replacement 

{^r});?^ - J{i, mtr!^ u^. ■ ■■u.y^ ■ ■ -^-N (158) 

in which we define 

= q.vu>'. (159) 

This approximation replaces each soft gluon entering the v4-jet by a fictitious 
gluon whose momentum and polarization are both in the w-direction. Before 
justifying the soft approximation, let us see what its consequences are. 

Once we make the soft approximation, each cut jet diagram is a contribu- 
tion to a product of matrix elements precisely of the form to which the Ward 
identity of Eq. (|133p can be applied, with the field $(x) now representing 
the field associated with the physically polarized parton which couples to the 
hard part. As a result, we have at out disposal a Ward identity, which can be 
used to factor soft lines from the jet subdiagram, by an iterative argument 
very similar to the one just used to factor longitudinally polarized collinear 
gluons from the hard part. The details of the argument are slightly more 
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a b 



Figure 26: Factorization of soft lines from a jet. 

complex because of the extra eikonal line, and we refer the interested reader 
to Refs. and |SFJ for details. Here we simply quote the result, which is 
illustrated in Fig. |2S1 and may be expressed as 

= ^ Jm{o,,' ^(^^ ^itc^'^'^£*{v^. rRi^^""^'^ (160) 

where Sac{v^, q1) stands for the lightlike eikonal line in the direction, 
to which have been connected those soft gluons to the left of cut C, with 
momenta gg and polarization indices cr(L). £* is defined similarly in terms 
of soft gluons to the right of the cut, with momenta and polarization 
indices a{R). Finally, we have made color indices a,b, ■ ■ ■ explicit, and d{R) 
is the dimension of the color representation of the physical parton. 
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The complete Green function may now be written in the form 



/ f ^ n / 1^ n / H ^'(^^f + 

xi^„,(g^e«■g^^^{g■})'"'^"^^ (lei) 

Notice that the jet function has now been factored from the rest of the 
process. When we now perform a sum over cuts, we can sum independently 
over the cuts of J and over the cuts of the rest of the diagram. In the rest of 
the diagram, we have a hard interaction producing the eikonal line £ and the 
jets in H. These are coupled by final state interactions with the soft gluons in 
5*. By the general reasoning of Sect. 17.61 these final state interactions cancel. 
Thus any leading region with soft exchanges cancels, and the factorization 
reasoning reverts to the arguments which apply in the scalar theory. We shall 
skip giving these details, and will close this subsection with a justification of 
the all- important soft approximation, Eq. ()158|) . 

The soft approximation consists of an approximation for the polarizations, 
and an approximation for soft momenta. The former may be justified by 
detailed power counting arguments |2S], but the underlying motivation is 
simply that gluon polarizations proportional to can couple to the A-jet 
by contracting into vectors proportional to p^, the momentum of the hadron 
A. Since u ■ p = p"^ is a large invariant, it will dominate by a power over 
invariants formed from the other internal momentum components present in 
J. Note, by the way, that in Feynman gauge gluon polarizations will have 
nonvanishing projections onto only if the soft subdiagram couples to other 
parts of the diagram as well as J. 

The approximation associated with the gluon momentum is more subtle. 
Recall that the A-jet is in the plus direction. We claim that one can neglect 
the transverse momentum of the gluon compared to its minus momentum. 
As we have seen in Sect. 18.41 this is nontrivial. In fact, regions where the 
transverse momentum is nonnegligible are leading by power counting. 

Recalling the one- loop discussion of Sect. 18.41 a typical denominator from 
the v4-jet on the left of the cut is of the form of 

(£ - Qif + ie^e^- 2tq~ + 21^ ■ q^, - Ig^J" + le (162) 

with ^^ a typical line momentum in the A-jet. We would like to set q^^ to 
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zero in all denominators like (|162j) . and all we need for that is 



(li I > 



■ qu - \qu 



(163) 



The relevant question is thus whether the q~ momentum contours are trapped 
at q~ = 0, at the scale of {2i±_ ■ — Note that poles of this type 

can only come from denominators from the v4-jet, and not from the hard part 
or the soft subdiagram. Now, although every jet line through which flows 
gives a pole at a position like ()162|) . close to the origin in the q~ plane, all of 
these poles are on the same side of the real axis. To see this, consider how 
each soft momentum q^ flows from the vertex where it attaches to J to the 
parton line that attaches J to H. In general, the q~ pole from any jet line is 
in the upper half plane if q^ flows in the opposite sense relative to the large 
plus momentum carried by that line. But may always be chosen to flow 
so that q~ is directed in this sense for each jet line on which it appears. This 
is evident from Fig. |^ Soft gluons to the right of the cut may be treated 
analogously. As a result, the q~ contours may all be deformed away from jet 
poles into a region where g^j may be neglected, and the soft approximation is 
justified along the deformed contour. By Cauchy's theorem, it is also justified 
in the original integral. Thus, the factorization program may be carried out 
in e~^e~ annihilation. 

9.3 The Drell-Yan cross section 

The thorniest factorization theorems involve two hadrons in the initial state. 
The Drell-Yan cross section for the process 



is the simplest of these, and has therefore received essentially all the atten- 
tion. Q'^ will represent the momentum of the lepton pair . The step 
from Drell-Yan to more complex processes, involving observed hadrons or 
jets in the final state is relatively straightforward, as indicated in Sect. 17.61 
Factorization for the Drell-Yan cross section has, at times, been the sub- 
ject of controversy [221 El 1121; although more recent work has, we believe, 
established its validity at all orders jSHESlEni- Nevertheless, as we shall ob- 
serve below, there is plenty of room for improvement in our understanding. 




(164) 
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The general leading region for the Drell-Yan process is shown in Fig. 
|22lb). After the sum over final states, all nonforward hadron jets are ab- 
sorbed into the hard subdiagram if, in the same way as in deeply inelastic 
scattering. In common with the deeply inelastic scattering and one particle 
inclusive e~^e~ annihilation processes, we can factor coUinear gluons from the 
hard part. Once this is done, the sum of cut Feynman diagrams for the 
Drell-Yan cross section is very similar to Eq. (jl57p for e^e~ annihilation, 

xH{Q^, U{u ■ Pa)v^, ■ pB)u''r' , (165) 

where the lightlike vectors = g'i^iU^ = have been chosen in the p'X,Pb 
directions, respectively and the parton momentum fractions are defined by 
= kA- u/pa ■ u and = ks ■ v/pb ■ v. Here Ja{U, Ui}) and Jb{Cb, {Qi}) 
are similar to the parton distribution function J defined for deeply inelastic 
scattering in Eq. (jlSOj) . except that soft gluons are still attached to them. 
Connections between the parton distribution and a soft subdiagram were 
absent in the deeply inelastic scattering cross section, because, after the 
sum over cuts, there was only one jet, which cannot by itself produce large 
invariants in numerator factors. In ()lfj5j) we have a soft subdiagram as in 
e+e" annihilation, but now interacting with the jets associated with the two 
incoming hadrons. 

Our basic problem is the same as in e^e~ annihilation, to show that con- 
tributions from any leading region with a nontrivial soft subdiagram cancel 
in the sum over final states and gauge invariant sets of diagrams. Then the 
remaining leading regions of Eq. ()165p are just of a form similar to Eq. ()144|) . 
and the arguments for factorization may be given as above for deeply inelas- 
tic scattering, eqs. ()152|) to p55|) . Naturally, we would like to proceed by 
analogy to e~^e~ annihilation. Thus, we would like to apply the soft approxi- 
mation to the jets, and factor the soft gluons from them. The jets would then 
contribute to parton distributions as in Eq. ()153p . and, once the remaining 
soft contributions cancel, we would derive the desired factorized form, Eq. 

The main obstacle to this program is shown in Fig. |27l which illustrates 
a typical low-order example. It shows a single soft gluon, q^, attached to the 
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Figure 27: Example illustrating obstacles to the soft approximation. 

A-jet. The soft momentum flows through two lines in the A-jet, an "active" 
jet line C. + q that carries positive plus-momentum into the hard part, and a 
"spectator" line p — £ — q that carries positive plus momentum into the final 
state. We saw in e+e^ annihilation that the criterion for the applicability of 
the soft approximation for qf^ is given by Ig"-^"*"! 3> |g_L ■ {2i± + q±)\, where 
£^ is any line in the jet along which g'^ may flow. But this condition may 
be satisfied along the entire q~ contour only if the contour is not pinched by 
poles on opposite sides of the real line. In e~^e~ annihilation they are not, 
but in Drell-Yan they are. This is illustrated by our example, since the poles 
in the q~ plane due to the jet propagators are approximately at 



which are on opposite sides of the contour, both at a distance of order 
2i± ■ q±/p'^ from the origin. Thus, in the Feynman integral associated with 
Fig. EZl the q~ contour is forced to go through a region in which the soft 
approximation fails, and we are unable to apply immediately the reasoning 
introduced for e^e~ annihilation. 

The resolution of this problem is rather technical, and may be found 
in Ref. (HH]- It may be understood most simply as a result of the Lorentz 
contraction of the colliding gluon fields, as in Sect. 18.11 In addition, we can 
give an intuitive picture here, based on semi-classical considerations in the 
center of mass frame. We consider the A-jet to be passing through a soft 
color field produced by the B-jet. Consider first the very softest part of the 
color field, with a spatial extent ~ (1 fm) GeV)]. The self-interactions 



{p - £f + 2(px - i^) ■ q^ - Ig^p + ie 

2{p - i) + 
-f + 2i^-q^ + \q^\^ -le 
2i+ 



(166) 
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of the partons in the A-jet are time dilated, but this field extends so far in 
space that it interacts with the partons on the same time scale as that of 
their self-interactions. However, on this distance scale, the soft gluons cannot 
resolve the hadron jet into individual partons. Thus the jet appears as a color 
singlet until the time of the hard interaction, at which time it acquires color 
because one parton is annihilated. The result is that the only interactions of 
the very soft color field occur long after the hard scattering event, and such 
interactions cancel because of unitarity. Consider now the part of the color 
field of smaller spatial extent, say 1 fm. The point is that the interactions of 
this color field with the spectator partons in the A-jet don't really matter. 
The reason is that the self-interactions of the partons in the A-jet are time 
dilated, so that the spectator partons do not interact with the active parton 
on the 1 fm time scale in which they interact with the color field. Since the 
spectator partons are not observed, unitarity implies that their interactions 
with the color field will not affect the cross section. As a technical trick, we 
could as well replace all of the spectator partons with an equivalent color 
charge located at = 0, right on top of the active quark. Then the color 
field sees a net color singlet in the initial state and a net colored charge in 
the final state. Again, the only interactions of the color field occur long after 
the hard scattering event, and such interactions cancel because of unitarity. 

Finally, note that the arguments given above are asymmetric between the 
two incoming jets. This is natural, because it is only necessary that one of 
the two incoming particles move at the speed of light for our arguments to 
apply. Indeed, factorization should hold in the (hypothetical) scattering of a 
truly lightlike particle with a massive particle at any center of mass energy. 
We should note that explicit two-loop calculations which show that infrared 
divergences cancel at leading power (although not at higher twist), have been 
carried out for the most part with one massive and one massless (eikonal) 
line 113 iU. 

10 Outlook and Conclusion 

In the foregoing, we have described the systematics of factorization for hard 
inclusive cross sections in QCD, and have discussed in some detail the nature 
of factorization proofs, first for (0^)6, and then for gauge theories. Along the 
way, we outlined a systematic approach to perturbative processes at high 
energy, based on the classification of leading regions. 
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As we have indicated above, the proof of factorization theorems in gauge 
theories is by no means a closed subject. Factorization proofs for inclusive 
processes is the first item of a whole list of subjects in which progress has been 
made, but for which important work remains to be done. In the following, 
we briefly discuss a few other significant topics which relate closely to the 
methods discussed in this chapter. Of great importance are extensions of the 
theorems to more general situations. 

10.1 Factorization Proofs 

Factorization proofs in nonabelian gauge theories have reached a certain level 
of sophistication in Refs. jSlI, [SHI and jHEl- Comparison with the discussion 
for however, shows that there is as yet in the literature no complete and 

systematic subtraction procedure in QCD of the type explained in Sect. 
even in the case of deeply inelastic scattering. A subtraction algorithm would 
eliminate any lingering uncertainty associated with overlaps between leading 
regions. Perhaps even more importantly, such a procedure should make it 
possible to develop bounds on corrections to leading power factorization the- 
orems, and to prove factorization theorems for nonleading power corrections, 
so-called "higher twist" . A model for this program is presumably to be found 
in the BPHZ formalism for deeply inelastic scattering cross sections in scalar 
and abelian gauge theories developed by Zimmermann suitably modi- 
fied to treat the extra infrared problems and gauge structure of QCD (see 
Sect. EI). 

It should also be noted that the Monte-Carlo event generators jlS] that are 
so widespread in analyzing data depend on generalizations of the factorization 
theorems; these generalizations have not yet gone significantly beyond the 
level of leading logarithms. 

In addition, we should mention that additional factorization theorems, of 
different but related forms, are central to the analysis of the elastic scattering 
of hadrons, which decrease as powers of the energy |46j . 

10.2 Factorization at Higher Twist 

It has been proposed |17j that generalized factorization theorems hold beyond 
the leading twist for a wide variety of cross sections. Most work on this 
possibility has been carried out for deeply inelastic scattering, where the 
systematics are best understood as a generalization of the operator product 
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expansion |48t H^l ED]- In particular, it has been shown that multiparton 
distributions may be defined in a natural way to parameterize soft physics 
at higher twist [13 EI] • 

In hadron-hadron scattering, factorization at higher twist is complicated 
by the infrared structure of perturbative QCD. We have seen in Sect. IHlthat 
leading twist factorization requires the cancellation of infrared divergences. 
It has been shown by explicit calculation > however, that infrared di- 

vergences do not cancel beyond a single loop in hadron-hadron scattering for 
QCD at higher twist. This is a sharp contrast between the nonabelian and 
abelian theories. At two loops, noncancelling divergences occur at the level 
m^/s^ in the Drell-Yan process. Refering to Sect. 18.11 this is precisely the 
level suggested by the classical relativistic kinematics of gauge fields. How 
one should interpret this lack of cancellation is not quite clear to us. The 
actual situation, including nonperturbative effects, may be better or worse 
than suggested by perturbative calculations The fact that perturbation 
theory respects factorization at m^/ s, however, makes it possible that factor- 
ization theorems may hold at this level, even for hadron-hadron scattering 

10.3 Factorization at the Boundaries of Phase Space 

A rich class of perturbative predictions involve the summation of corrections 
near boundaries of phase space in different processes. Near some of these 
boundaries, notably small Q± and small x, cross sections increase greatly. 

Along these lines, perhaps the most attention has been given to the 
Drell Yan cross section at measured transverse momentum dcr/dQ^d^Q^ 
inHEHlESl, with -C Q and the related two-particle inclusive cross section 
for e~^ + e~ — >■ A + B + X at measured transverse momentum j^EZlEHl- The 
complete leading-twist analysis of these cross sections begins with factorized 
forms of the type of Eq. ()165|) , in which soft partons have been factored from 
jets, but not yet cancelled. At the boundary of phase space the cancellation 
of soft gluons outlined in Sect. El still occurs, but is incomplete. All infrared 
divergences still cancel at leading power, but finite remainders depend on the 
small parameter in the problem, for instance the transverse momentum in 
the cross sections cited above. By developing generalizations of the renor- 
malization group equation for each of the functions in the factorized form 
Eq. (jl65p |4Uj . it is possible to resum systematically higher order corrections 
to these quantities. 
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This general approach can be apphed in a number of other physically 
important situations. For instance, the r = Q^/s 1 limit in the inclusive 
Drell-Yan cross section is related to the normalization of the Drell-Yan cross 
section, the "K-factor" jSHl 1101 lEH 102] • It is possible to sum corrections 
which are singular at r = 1 fT^ lU^ . An interesting feature of the result 
is that it is sensitive to high orders in perturbation theory [01] through the 
running coupling. Because of this, it gives a measure of the sensitivity to 
higher- twist effects of perturbative predictions based on factorization jOO]- 
This sensitivity is found to be nonnegligible in some, but not all, regions of 
physical interest. 

Another regime, which is of crucial importance for experiments at the 
Tevatron and SSC, is the x — limit in hadron-hadron scattering. The cross 
sections get into the range of tens of millibarns, which is enormous compared 
to typical cross sections at larger x. So far, much work has concentrated 
on the behavior of parton distributions at small x [001 1001 lOZl lOH] , assuming 
the vahdity of the standard factorization formulas Q and (fTTj) . From a 
more general point of view, factorization has been shown to hold explicitly 
in leading logarithms in x jOH]- We would like to suggest, however, that 
factorization theorems need a more extensive examination in this region. 

In conclusion, we emphasize that essentially every calculation in pertur- 
bative QCD is based on one factorization theorem or another. In view of 
this, progress toward developing perturbative QCD as a quantitative system 
requires further understanding of the systematics of factorization. 
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